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1.0 


INTRODUCTION 


In  the  last  25  years,  a  great  deal  of  progress  has  been  made 
in  the  design  of  the  helicopter.  Rotary  wing  aircraft  are  now 
faster,  more  efficient,  quieter  and  carry  more  payload  than  ever 
before.  For  example,  in  just  one  generation  hover  efficiency 
(figure  of  merit)  has  been  increased  by  about  15%,  which  trans¬ 
lates  into  approximately  double  the  payload  capability.  These 
advances  are  a  direct  result  of  improved  understanding  of  heli¬ 
copter  aerodynamics  that  has  come  with  the  development  of  better 
analytical  methods  and  experimental  techniques.  The  theoretical 
methods  available  today  are  certainly  much  more  advanced  than 
those  available  in  the  1950 *s. 


However,  it  must  be  surprising  for  those  of  us  who  have  been 
involved  in  the  development  of  methods  for  the  prediction  of 
rotary  wing  airloads  to  realize  that  it  was  approximately  25 
years  ago  that  Gray!'2  first  characterized  the  hovering  rotor 
wake  as  composed  of  strong,  rolled  up  helical  tip  vortices  and 
separate  inboard  sheets  of  weak  vorticity.  Since  that  time, 
research  and  development  efforts  have  concentrated  on  the  pre¬ 
scribed  wake  method.  Still,  it  was  almost  15  years  before 
Landgrebe2  had  parameterized  the  prescribed  wake  constants  in 
terms  of  blade  geometric  properties  and  rotor  thrust  coefficient 
based  on  a  series  of  model  rotor  tests  and  a  lifting-line 
representation  of  the  rotor  blades.  During  the  1970's,  additional 
model  rotor  tests  have  led  to  continued  refinements  to  the  pre¬ 
scribed  wake  constants,  and  maximum  blade  circulation  has  replaced 


2 
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■^Gray,  R.B.,  "On  the  Motion  of  the  Helical  Vortex  Shed  from 
a  Single-Bladed  Hovering  Helicopter  Rotor  and  its  Application  to 
the  Calculation  of  the  Spanwise  Aerodynamic  Loading",  Princeton 
University  Aero.  Engr.  Dept.,  Report  No.  313,  September  1955. 

2Gray,  R.B.,  "An  Aerodynamic  Analysis  of  a  Single-Bladed 
Rotor  in  Hovering  and  Low-Speed  Forward  Flight  as  Determined 
from  Smoke  Studies  of  the  Vorticity  Distribution  in  the  Wake", 
Princeton  University  Aero.  Engr.  Dept.,  Report  No.  356,  September 
1956. 


^Landgrebe,  A. J. ,  "An  Analytical  and  Experimental  Investiga¬ 
tion  of  Helicopter  Rotor  Hover  Performance  and  Wake  Geometry 
Characteristics" ,  USAAMRDL  Tech.  Rept.  71-24,  Eustis  Directorate, 
U.S.  Army  Air  Mobility  Research  and  Development  Laboratory,  Fort 
Eustis,  Virginia,  June  1971,  AD  728835. 
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rotor  thrust  in  the  empirical  prescript ion . 4 ' ^  That  these 
analytical  methods  demonstrate  good  correlation  for  integrated 
loads  for  a  large  number  of  conventional  rotors  should  be 
expected  since  the  prescribed  wake  constants  are  intimately 
related  to  the  theoretical  methods  themselves.  However,  these 
methods  are  less  successful  when  compared  with  known  collective 
settings, ^  which  implies  that  although  integrated  performance  is 
predicted  correctly,  local  sectional  loads  are  probably  not 
properly  calculated.  Certainly,  if  the  wake  itself  deviates 
from  the  experimental  data  base  for  the  prescribed  wake  constants, 
performance  predictions  could  be  significantly  in  error.  In 
fact,  detailed  performance  predictions  for  the  new  technology 
rotor  blade  tips  that  exist  today  still  remain  uncertain. 


The  lifting-surface  method  reported  in  Reference  7  should 
help  in  analyzing  the  new  plan  forms,  providing  the  linearized 
lifting-surface  assumptions  are  met  and  the  prescribed  contrac¬ 
tion  boundaries  are  acceptable.  The  resulting  computer  program, 
HOVER,  is  described  in  detail  in  Reference  8. 


Landgrebe,  A.J.  ,  Moffit,  R.C.,and  Clark,  D.R.,  "Aerodynamic 
Technology  for  Advanced  Rotorcra f t — Part  I",  J.  Am.  Hel .  Soc., 
Vol.  22,  No.  2,  April  1977. 

^Kocurek,  J.D.,  "Hover  Performance  Methodology  at  Bell 
Helicopter  Textron",  36th  Annual  Forum  of  the  American  Helicopter 
Society,  Preprint  No.  80-3,  May  1980. 

^Caradonna,  F.X.,and  Tung,  C.,  "Experimental  and  Analytical 
Studies  of  a  Model  Helicopter  Rotor  in  Hover",  Sixth  European 
Rotorcra ft  and  Powered  Lift  Aircraft  Forum,  Paper  No.  25, 

Bristol,  England,  September  1980. 

^Summa,  J.M.,and  Clark,  D.R.,  "A  Lifting-Surface  Method  for 
Hover/Climb  Airloads",  35th  Annual  National  Forum  of  the  American 
Helicopter  Society,  Preprint  No.  79-3,  May  1979. 

^Summa,  J.M.,  "Evaluation  of  Blade  Tip  Planform  Effects  on 
Hover  Performance" ,  Contract  DAAJ02-7 6-C-0069 ,  Applied  Technology 
Laboratory,  U.S.  Army  Research  and  Technology  Laboratories 
(AVRADCOM) ,  Fort  Eustis,  Virginia,  to  be  published. 


With  this  program,  the  axial  displacements  of  the  tip  vor¬ 
tices  are  relaxed  to  force-free  locations  while  the  contraction 
boundary  is  prescribed  as  in  the  earlier  methods. Convergence 
of  the  relaxed  wake  and  rotor  loading  with  the  program  have  been 
demonstrated  for  a  small  number  of  rotors,  but,  unfortunately, 
extensive  correlation  efforts  and  calculations  for  newer  blade 
geometries  have  not  been  undertaken.  It  is  anticipated  that  such 
calculation  will  be  completed  within  the  next  year. 

Meanwhile,  Green's  function  surface  singularity  method, 
described  in  this  report,  has  been  developed  to  examine  the  de¬ 
tailed  pressure  distributions  on  rotor  blades  of  arbitrary 
shape,  and  to  include  the  thick  blade  surface  vorticity  effects  at 
first  blade  passage.  Of  course,  this  new  method  is  also  the  first 
step  toward  the  eventual  goal  of  the  elimination  of  the  two-dimen¬ 
sional  constraints  in  the  profile  power  predictions. 

The  basic  methodology  for  the  surface  singularity  code  was 
developed  for  fixed-wing  applications  involving  unsteady  oscilla¬ 
tions  and  has  been  described  in  References  9  and  10.  Rotary  wing 
wake  technology  has  been  added  to  the  basic  method  along  with  the 
necessary  methods  to  calculate  rotorcraft  airloads.  The  result¬ 
ing  computer  code,  named  ROTAIR,  is,  therefore,  based  on  tech¬ 
nologies  developed  in  HOVER  and  the  fixed-wing  computer  programs. 

This  report  describes  briefly  the  basic  fundamentals  of  the 
surface  potential  method  as  applied  to  the  rotor  in  hover  or 
climb.  Further  details  of  the  fundamentals  can  be  found  in 
References  9  and  10.  The  complete  description  of  the  blade  sur¬ 
face  geometric  description  in  ROTAIR  is  described  in  Section  3.0, 
while  the  wake  modelling  and  geometry  are  presented  in  Section 
4.0.  Rather  than  repeating  technology  developments  already  re¬ 
ported  in  great  detail  in  Reference  8,  wake  modi f ications  made 
to  the  methods  for  incorporation  into  ROTAIR  are  emphasized. 
Finally,  the  method  of  rotor  airloads  computation  and  data  cal¬ 
culations  are  discussed  in  Sections  5.0  and  6.0,  respectively. 


9 

Maskew,  B. ,  "Unsteady  Potential  Flow  Analysis  of  Rotor 
Blade  Tip  Shapes",  Interim  Report  and  Program  User  Guide,  Contract 
NAS1-15472,  NASA  Langley  Research  Center,  April  1980. 

l^Maskew,  B. ,  "Influence  of  Rotor  Blade  Tip  Shape  on  Tip 
Vortex  Shedding--An  Unsteady,  Inviscid  Analysis",  Presented  at 
the  36th  Annual  Forum  of  the  American  Helicopter  Society, 

Preprint  No.  80-6,  May  1980. 
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2.0  SURFACE  SINGULARITY  METHOD  FOR  ROTORS 


The  complete  mathematical  formulation  for  the  general  motion 
of  multiple  bodies  through  an  ideal  fluid  was  rigorously  developed 
in  Appendix  A  of  Reference  8.  In  that  formulation,  the  boundary 
condition  of  zero  flow  penetration  at  the  blade  surface  was  used 
to  obtain  solutions  for  the  blade  surface  doublet  distribution. 
Consequently,  the  Neumann  velocity  boundary  condition  of  no  flow 
across  the  boundary  was  enforced  on  the  external  (outside  of  the 
body)  flow  field,  and  a  doublet  only  distribution  was  utilized. 
However,  for  the  motion  of  closed  bodies,  the  internal  Dirichlet 
boundary  condition  of  zero  potential  inside  the  body  can  be  used 
equally  as  well.  In  fact,  there  are  numerical  advantages  to 
solving  for  the  velocity  potential  directly,  since  it  is  one 
order  less  singular  than  the  velocity  induced  by  a  potential 
doublet  distribution.  For  example,  in  Reference  11  it  was  demon¬ 
strated  that  for  comparable  density  of  control  points  where  the 
boundary  conditions  are  satisfied,  the  low-order  method  gives 
comparable  accuracy  to  the  higher-order  solutions.  It  was  also 
shown  that  problems  associated  with  some  earlier  low-order  panel 
methods,  e.g.,  leakage  in  internal  flows  and  junctions  and  also 
poor  trai ling-edge  solutions,  do  not  appear  for  the  velocity 
potential  method.  Consequently ,  the  surface  singularity  method 
for  the  rotor  problem  uses  the  Dirichlet  boundary  condition  for 
the  velocity  potential.  For  completeness,  the  mathematical  for¬ 
mulation  given  in  Reference  9  is  repeated  here  with  only  slight 
modi f icat ions  for  the  rotor  problem  of  concern  here.  Also,  the 
formulation  for  a  single  blade  only  is  explicitly  developed  in 
the  following  sections  since  the  influences  of  secondary  blades 
and  wakes  are  obtained  by  simple  rotations  for  the  case  of  hover 
or  climb.  The  reader  is  referred  to  Reference  7  for  details  of 
the  multiple-blade  problem. 


11-Maskew,  B.,  "Prediction  of  Subsonic  Aerodynamic  Character- 
istics--A  Case  for  Low-Order  Panel  Methods",  AIAA  19th  Aerospace 
Sciences  Meeting,  AIAA-8 1-0252 ,  January  1981. 


2.1 


BASIC  FORMULATION 


Consider  the  whole  of  space  divided  into  two  regions  by 
the  surface  of  the  blade  and  assume  the  existence  of  Laplacian 
velocity  potential  distributions  in  the  two  regions,  i.e.,  *  in 
the  flow  field  and  in  the  blade  interior.  If  we  now  apply 
Green* s  third  identity  to  the  two  regions,  then  the  total  poten¬ 
tial  at  a  point,  P,  on  the  inside  surface  of  the  blade  can  be 
written: 


4"rp  '//(*-♦  i)s  ■  '•’(?)  OS  -  2,(»  - 

s-p  p 


+  //  (  fu  -  "l)-  •  •  (?)  ow 
w 


-//  ?  2  •  (v:  -  7‘i)ds  +  4’f-r  (1» 

s 


Here,  r  is  the  length  of  the  vector  from  the  surface  element  to 
the  point,  P,  and  S-P  signifies  that  rhe  point,  P,  is  excluded 
from  the  surface  integration.  Equation  (1)  includes  the  contri¬ 
bution  from  the  wake  surface,  W. 

The  Dirichlet  boundary  condition  is  now  applied  in  the 
interior  region  to  render  a  unique  distribution.  Two  forms  have 
been  considered  in  the  present  work:  (1)  J>i  =  0  and  (2)  = 

The  first  form  gives  a  simpler  and  more  cost-effective  method 
for  the  qeneral  case,  but  the  second  form,  which  is  used  for  the 
present  calculations ,  has  proven  to  be  more  reliable  in  practice 
For  the  second  form,  Eq  (1)  becomes 
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n  *  jdS  (2) 
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where  ^  ,  the  perturbs 1  i on  potential  in  the  flow  field,  has  been 
substituted  for  v  -  $  . 

iV> 

The  first  two  terms  in  Eq  (2)  give  the  perturbation 
potential  due  to  a  distribution  of  normal  doublets  of  strength, 
on  the  blade  surface.  Similarly,  the  third  term  represents  a 
doublet  distribution  of  strength,  <J>u  -  <J>L,  on  the  wake  and  the 
fourth  term  represents  a  source  distribution  of  strength,  n  •  V<j>, 
on  the  blade  surface. 

Equation  (2)  is  basically  the  same  as  the  formulation  given 
by  Morino,  12  who  used  a  direct  application  of  Green's  theorem 
in  the  flow  field.  The  present  approach  to  the  problem  is  a 
special  case  of  a  multi-domain  formulation  which  has  led  to  a 
more  general  three-dimensional  method  in  which  large  regions  of 
separated  flow  are  modelled  in  a  similar  way  to  that  in  the  CLMAX 
program.1-3 

The  source  term  in  Eq  (2)  can  be  evaluated  directly  from 
the  condition  of  no  flow  penetration  at  the  surface.  The  flow 
velocity  relative  to  the  blade-fixed  frame  is 

Y  =  v  -  V  -  w  .  R,  (3) 


12 

Morino,  L. ,  Chen,  L.T.,and  Socio,  E.O.,  "Steady  and  Oscil¬ 
latory  Subsonic  and  Supersonic  Aerodynamics  Around  Complex  Con¬ 
figurations",  AIAA  J.  ,  Vol.  13,  No.  3,  March  1975,  pp.  368-374. 

13 

Maskew,  B.,  and  Dvorak,  F.A. ,  "The  Prediction  of  C 
Using  a  Separated  Flow  Model",  J.  Amer.  Hel.  Soc.  ,  April  1978. 


where  the  perturbation  velocity,  v  =  -V$. 

For  zero  penetration,  V  •  n  =  0.  Hence, 


n  •  =  -n  •  -  n  •  a)  R, 


and  Eq  (2)  becomes 


(4) 


This  is  the  basic  equation  of  the  method.  It  is  solved  for 
the  unknown  surface  perturbation  potential,  t ,  or  surface  doublet 
distribution  for  a  given  wake  geometry.  The  doublet  distribution, 

-  (^t,  on  the  wake  is  equal  to  the  surface  doublet  jump  at  the 
blade  trailing  edge,  as  specified  by  the  Kutta  condition.  The 
wake  geometry  is  either  prescribed  or  relaxed  to  a  force- free  loca¬ 
tion  by  the  methods  described  in  Reference  8.  Once  the  surface 
potential  distribution  is  known,  the  surface  velocities  are  ob¬ 
tained  by  differentiation,  and  the  blade  surface  pressures  are 
calculated  in  the  usual  manner. 


2 . 2  NUMERICAL  PROCEDURE 


The  general  arrangement  of  the  configuration  is  shown 
in  Figure  1.  The  x,y,z  coordinate  system  with  unit  vectors,  _i, 
j^,k,  is  fixed  relative  to  the  blade.  For  the  present  applications, 
the  blade  rotates  about  the  z-axis  at  a  rate,  Q ,  and  climbs  at  a 
speed,  Wc ,  along  the  z-axis. 
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EDGE  SEPARATION  LINE 
PRESCRIBED 


Figure  1.  General  Arrangement  of  the  Configuration. 


The  numerical  procedure  has  options  for  both  a  prescribed 
and  relaxed  wake  calculation  scheme.  The  flow  chart  for  the 
procedure  is  shown  in  Figure  2.  The  input  data  is  assembled  in 
three  main  parts:  basic  data,  surface  patch  geometry  description , 
and  basic  wake  geometry  description.  Details  of  the  input  param¬ 
eters  and  the  program  output  are  given  in  Appendixes  B  and  C  as 
well  as  a  more  detailed  flow  diagram  showing  the  connections  be¬ 
tween  subroutines  in  the  code. 

The  surface  of  the  blade  is  represented  by  planar  quadrila¬ 
teral  panels  over  each  of  which  the  doublet  and  source  distribu¬ 
tions  are  assumed  constant.  Special  features  are  included  in  the 
present  code  to  extend  the  panelling  around  the  tip  edge  and  to 
allow  the  user  to  extend  the  wake  along  the  tip  edge  (Figure  1) . 
Automatic  panelling  routines  are  installed  to  simplify  user  input. 
The  geometry  description  is  given  in  detail  in  Section  3.0. 

Equation  (4)  is  satisfied  simultaneously  at  a  point  at  the 
center  of  each  panel.  If  there  are  N  panels  representing  the 
blade  surface,  Eq  (4)  becomes 


where  %  is  the  number  of  panels  in  the  wake,  and  =  Wck. 
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Figure  2.  Flow  Diagram  for  the  Method. 
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is  the  source  distribution  due  to  blade  rotation  about  th 
pitch  axis,  and 


I  N  I 

!  £  2K  Bjk  1  '« 


are  the  .  o^ronents  of  a  three-part  source  distribution  due  to 
the  ^dative  translation  of  the  blade  and  the  onset  flow.  (Note: 
In  the  present  case,  only  the  z-component  is  nonzero  for  the 
climb  condition.) 

The  quantities  and  Cjk  are  the  velocitv  potential  in¬ 

fluence  coefficients  for  the  constant  source  and  doublet  distri¬ 
butions,  respectively,  on  panel  K  acting  on  the  control  point  on 
panel  J.  These  include  contributions  from  the  respective  second¬ 
ary  blade  panels.  Expressions  for  these  influence  coefficients 
have  been  given  by  Morino  in  Reference  12  based  on  hyperbolic 
paraboloidal  panels.  Slightly  different  expressions  are  installed 
in  the  present  code  based  on  planar  panels. 

Equation  (5)  is  solved  by  a  direct  method  for  N  £  320  and  by 
an  iterative  method  for  N  >  320. 

The  surface  gradient  of  M  is  evaluated  on  each  panel  by 
dif ferentiating  a  two-way  parabolic  fit  through  the  doublet  values 
on  the  panel  and  its  four  immediate  neighbors.  At  the  separa¬ 
tion  lines  a  simple  differencing  is  applied  for  the  gradients 
approaching  the  separation  line.  The  surface  pressure  distribu¬ 
tion  and  resultant  airloads  are  calculated  as  described  in  Sec¬ 
tion  5.0. 

The  prescribed  wake  iteration  loop  requires  that  the  user 
input  the  empirical  wake  constants  for  the  desired  rotor  thrust 
coefficient  and  an  approximate  radial  position  of  the  maximum 
circulation.  The  iteration  continues  until  the  maximum  circula¬ 
tion  position  is  unchanged.  Once  the  collective  iteration  loop 
is  added  to  the  program,  the  blade  collective,  075/  would  then 
be  updated  until  the  calculated  thrust  equals  the  requested 
thrust.  If  requested,  a  relaxed  wake  iteration  loop  similar  to 
that  described  in  Reference  8  is  then  calculated  for  a  fixed- 
blade  collective. 
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3.0  BLADE  GEOMETRY  DEFINITION 


This  section  describes  the  way  the  geometry  routines  pro¬ 
cess  user  input  to  generate  the  discrete  panel  representation 
of  the  rotor  blade.  Details  are  included  to  emphasize  the 
versatility  of  the  routines.  With  the  exception  of  the  method 
of  providing  for  coning  angle  rotation,  8,  and  collective  pitch 
rotation,  675,  the  geometry  description  is  carried  over  from 
Reference  9. 


3.1  BLADE  SURFACE  ORIENTATION 


When  defining  the  surface  geometry  of  a  rotor,  the  primary 
blade  is  described  in  its  own  local  coordinate  system  for  con¬ 
venience.  The  blade  surface  is  then  relocated  to  the  appropriate 
position  with  respect  to  the  rotor  hub. 

The  blade  specification  starts  with  the  appropriate  trans¬ 
formation  information  which  converts  from  the  blade  definition 
system,  or  "blade  coordinate  system"  (B.C.S.)  to  the  hub  coordi¬ 
nate  system  (H.C.S.).  This  information  includes  the  flapping 
hinge  radial  position,  YFLAP ,  the  blade  coning  angle,  8/  and  the 
blade  collective  setting,  n75.  As  illustrated  in  Figure  3,  the 
oriqin  of  the  B.C.S.  is  positioned  at  YFLAP  distance  from  the  axis- 
o.f-rotation  on  the  yh.C.S.  axis ,  the  blade  (or  B.C.S.)  is  then 
coned  8  degrees  about  the  axis,  and  finally  rotated  to 

the  requested  ^75  degrees  about  the  Yp  c  S  axis.  Since  the  rotor 
radius  is  used  for  nondimensionalization’,  care  must  be  exercised 
to  ensure  that  the  final  blade  radius  includes  Y^LAP .  A  test  is 
made  at  the  geometry  stage  in  ROTAIR,  and  if  in  error,  the  pro¬ 
gram  is  halted  with  a  message  for  the  user  to  check  the  radius 
defining  stations.  Finally,  the  blade  surface  is  described  with¬ 
in  the  B.C.S.  in  patches  of  panels. 

3.2  PATCHES 

3.2.1  Convention 

Panels  representing  the  surface  of  the  configuration  may  be 
assembled  in  a  number  of  PATCHES.  Each  patch  has  a  regular  array 
of  panels  arranged  in  rows  and  columns  (Figure  4) .  A  patch  is 
basically  a  four-sided  shape  when  "opened  out".  It  must  always 
be  regarded  as  such,  even  if  one  of  the  sides  or  two  opposite 
sides  are  made  zero,  or  if  some  of  the  sides  have  kinks.  In  the 
following  discussions,  patches  will  often  be  regarded  as  rectangu¬ 
lar;  this  is  purely  a  convenience  for  discussing  relationships  and 
is  not  a  shape  restriction.  Our  view  of  the  patch  will  always  be 
from  the  outside;  i.e.  ,  looking  onto  the  surface  from  a  point  in 
the  flow  field. 
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Figure  4.  Patch  Conventions. 


For  convenience  the  terms  "chordwise"  and  "spanwise"  are 
used  to  describe  the  directions  of  the  panel  columns  and  rows, 
respectively  (Figure  4) .  These  directions  are  analogous  to  the 
conventional  wing  layout,  but  in  the  patch  context  these  direc¬ 
tions  are  not  restricted  to  the  x  and  y  directions,  respectively. 

For  example,  on  a  patch  representing  the  wing  tip  it  is  con¬ 
venient  to  have  the  columns  of  panels  vertical  and  the  rows  of 
panels  to  be  parallel  to  the  wing  chord;  therefore,  in  this  case 
the  "chordwise"  direction  on  the  patch  is  actually  vertical  while 
its  "spanwise"  direction  is  along  the  wing  chord  (see  Figure  10) . 

Patch  geometry  is  defined  using  chordwise  lines  called  SECTIONS. 
(These  are  described  later  in  3.2.2.)  A  set  of  sections  distri¬ 
buted  spanwise  across  a  patch  defines  the  patch  surface.  The  con¬ 
vention  adopted  here  is  that  points  defining  a  section  shape  pro¬ 
ceed  from  top  to  bottom  (Figure  5) .  (In  the  case  where  a  patch 
represents  the  main  surface  of  a  rotor  blade,  this  convention 
causes  the  points  defining  each  section  to  proceed  from  the 
trai ling-edge  lower  surface  and  finish  at  the  trai ling-edge 
upper  surface.)  In  our  view  of  the  patch,  the  order  of  the  sec¬ 
tions  proceeds  in  the  positive,  spanwise  direction  (Figure  5) . 

For  the  purpose  of  automatically  connecting  panels  from  one 
patch  to  another,  it  is  important  to  identify  patch  sides.  The 
convention  adopted  here  is  that  the  first  and  last  sections  de¬ 
fining  a  patch  correspond  to  sides  1  and  3,  respectively.  With 
this  convention,  the  order  of  the  sides  is  anticlockwise  (Figure 
4) .  The  order  of  the  corner  points  follows  the  same  sequence  as 
the  sides,  starting  with  1  at  the  top  of  side  1.  All  the  panels 
within  a  patch  take  the  same  side  and  corner  point  convention  as 
for  the  patch.  For  convenience,  the  panels  are  referred  to  in 
ROWS,  which  run  spanwise,  and  COLUMNS,  which  run  chordwise  (Figure 
4) .  Panel  arrangements  within  a  patch  are  referred  to  by  ROWS  x 
COLUMNS . 


3.2.2  Sections 

Each  section  of  a  patch  may  be  defined  in  its  own  local 
coordinate  system,  referred  to  as  the  section  coordinate  system, 
or  S.C.S.  The  user  provides  the  necessary  information  to  trans¬ 
form  the  S.C.S.  into  the  B.C.S.  at  the  beginning  of  each  section. 
This  transformation  is  performed  immediately  after  a  section’s 
geometric  description  is  complete.  This  transformation  is  sepa¬ 
rate  from  that  described  earlier  in  3.1  in  which  the  complete 
blade  qeometry  is  converted  into  the  H.C.S.  (at  which  stage  the 
S.C.S.  geometry  is  discarded).  This  double  transformation 
(both  1  evels  of  which  are  optional)  offers  useful  flexibility 
when  preparing  the  input  data.  One  particular  advantage  is  that 
the  geometric  relationships,  especially  the  rotations,  are  kept 
reasonably  simple  without  sacrificing  qeneralitv. 
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The  information  required  to  transform  from  the  S.C.S.  into 
the  B . C . S .  (see  Figure  6)  consists  of  :  (1)  the  translation 

vector  (STX,  STY,  STZ) ,  which  is  the  position  vector  of  the  S.C.S. 
origin  expressed  in  the  B.C.S.  coordinates;  (2)  a  scaling  factor 
(default  1.0),  which  is  applied  in  the  S.C.S.;  (3)  the  rotation 

anqle  (a,  degrees)  about  the  y-axis  of  the  S.C.S.;  and  (4)  the 
angle  (0,  degrees)  in  the  B.C.S.  x-y  plane,  between  the  projection 
of  the  S.C.S.  x-axis  and  the  B.C.S.  x-axis.  For  most  rotor  blades, 
the  9  rotation  will  be  zero  and  the  section  twist  anqle,  in  this 
case  ot,  will  be  used  to  specifv  aribtrary  radial  twist  distribu¬ 
tions  (Fiqure  7)  .  It  is  not  necessarv  for  the  resultant  blade 
twist  anqle  at  the  3/4  radius  to  be  zero.  The  qeometry  packaqe 
computes  the  basic  twist  of  the  defined  blade  surface  at  the  3/4 
radius  station  and  rotates  the  blade  to  the  specified  O75.  How¬ 
ever,  because  of  the  rotation  procedures  described  in  3.1  for  the 
final  blade  position,  the  translation  vector  (STX,  STY,  STZ)  of 
the  blade  root  section  must  be  set  so  that  the  Y]$mc.S.  axis  passes 
through  the  root  section  quarter  chord  point. 

The  contour  line  of  each  section  is  defined  by  a  set  of 
BASIC  POINTS  (BX,  BZ).  These  points  may  be  used  as  panel  corner 
points;  i.e.,  MANUAL  PANELLING,  in  which  case  the  user  must  take 
care  over  the  number  of  input  points.  Alternatively,  an  AUTO¬ 
MATIC  PANELLING  ROUTINE,  referred  to  as  the  A.P.R.,  may  be  acti¬ 
vated,  which  interpolates  through  the  basic  points  to  form  a  new 
set  of  points  corresponding  to  panel  corner  points.  (This  is  just 
a  temporary  set,  as  the  user  may  opt  to  use  the  A.P.R.  in  the 
spanwise  direction  as  well,  in  which  case  the  section  points  do 
not  necessarily  line  up  with  panel  edges.) 

Several  options  which  provide  great  f  1  * i  Pilit  j*  .Oen  pre¬ 
paring  the  input  have  been  provided  for  de'  '  :ung  the  basic  points. 
The  options  may  be  exercised  at  the  section  level  so  the  input 
form  mav  be  chanqed  from  section  to  section.  The  options  available 
at  this  time  are  described  below  and  are  controlled  by  the  value 
of  INMODE.  INMODE  values  of  1  throuqh  4  are  illustrated  in 
Fiqure  8. 

INMODE  values  of  1,  2  and  3  are  used  when  a  section  lies  in 
one  of  the  reference  planes  of  the  chosen  S.C.S.;  in  these  cases 
we  have  a  constant  coordinate,  x,  y,  or  z,  respectively.  With 
one  coordinate  fixed  we  need  input  only  two  coordinates  for  each 
basic  point,  e.g.,  y  and  z  when  INMODE  =  1.  Provision  is  made 
to  specify  a  third  quantity  to  give  a  local  adjustment  to  the 
"constant"  coordinate,  e.g.,  when  using  INPUT  =  2  we  may  specify 
x,  z  and  £y.  Usually  the  <5  quantity  is  left  blank  (i.e.,  0)  ; 
however,  nonzero  S  values  allow  the  user  to  "bend"  the  input 
section  about  the  specified  "constant"  coordinate.  The  basic 
value  for  the  constant  coordinate  is  zero  until  the  section  points 
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Figure  8.  Basic  Point  input  Options  1  through  4. 
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are  transformed  into  the  B.C.S.,  so  the  value  of  that  coordinate 
in  the  B.C.S.  must  be  provided  in  the  transformation  information. 
INMODE  value  of  4,  which  requires  all  three  components  of  each 
basic  points  position  vector,  is  used  when  defining  a  completely 
arbitrary  section  shape.  Finally,  if  the  section  shape  has  a 
NACA  4-digit  form,  then  the  INMODE  parameter  can  be  set  to  5 
and  the  code  generates  the  panelling  directly.  The  trailing- 
edge  thickness  has  been  made  zero  by  modifying  the  coefficient 
in  x4  in  the  4-digit  equation,  and  the  coordinates  are  generated 
in  the  INMODE  =  2  format ,  i.e.,  x,  z  with  y  =  0. 

Zero  or  negative  values  allow  the  present  section’s  basic 
points  to  be  copied  over  completely  from  any  previously  defined 
section.  The  section  number  is  (-INMODE)  except  when  INMODE  =  0; 
the  latter  copies  over  the  points  from  the  section  just  completed. 
The  section  number  specified  is  the  absolute  number  from  the 
beginning  of  the  input  and  includes  other  copied  sections  as  well 
as  sections  which  may  have  been  generated  automatically.  If  the 
section  countinq  becomes  complicated,  alternative  ways  of  copying 
are  available  as  described  later  in  3.2.5.  The  basic  points  are 
copied  from  the  S.C.S.  set  (i.e.,  as  originally  specified)  and  are 
then  transformed  to  the  present  B.C.S.  according  to  the  new 
section1 s  transformation  information. 


3.2.3  Chordwise  Regions 

The  basic  points  defining  a  section  may  be  assembled  in 
a  number  of  CHORDWISE  REGIONS  for  the  purpose  of  controlling 
the  panel  density  and  distribution  on  that  section.  In  addition, 
the  option  on  manual  or  automatic  panelling  is  selected  at  the 
chordwise  region  level,  allowing  the  user  to  switch  from  one  to 
another  within  each  section  wherever  he  chooses.  Chordwise 
regions  are  used  only  as  an  input  convenience  and  are  discarded 
in  the  programs  as  soon  as  the  surface  panelling  is  complete, 

A  chordwise  region  must  end  on  a  basic  point  called  a  NODE 
POINT(Figure  9) .  A  NODE  CARD,  containing  the  chordwise  region 
panelling  information  (see  below) ,  inserted  after  a  basic  point 
in  the  input  deck  identifies  that  point  as  the  end  of  a  chord- 
wise  region.  Node  points  are  usually  placed  at  "problem”  areas 
where  large  velocity  gradients  are  expected  to  occur,  e.g.,  flap 
hinge  line,  leading  edge,  close-interference  regions,  but  the 
user  can  place  them  wherever  he  wishes  to  change  from  one  panel 
scheme  to  another.  Four  types  of  node  points  are  provided  at 
this  time  and  are  described  below. 
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BASIC  POINTS  DEFINE 
SECTION  SHAPE 


SEQUENCE  OF 
BASIC  POINTS 
ON  EACH  SECTION 


Figure  9. 


Chordwise  Regions  on  a  Section. 


The  information  on  a  NODE  CARD  consists  of  just  three  integers: 

(i)  NODEC  identifies  the  node  point  and  its  type. 

(ii)  NPC  is  the  number  of  panels  to  be  generated  by  the 

A.P.R.  in  the  chordwise  region  just  completed — 
a  zero  value  gives  manual  panelling. 

(iii)  TNTC  controls  the  form  of  the  distribution  in  the 
automatic  oanellina  mode  and  is  inactive  in 
the  manual  oanellinq  mode. 

(The  C  on  the  end  of  each  quantity  distinguishes  the  chordwise 
from  the  corresponding  spanwise  quantities,  which  end  in  S; 
see  Section  3.2.4) 

NODEC  values  of  1  or  2  specify  the  end  of  a  chordwise  region 
with,  respectively,  continuous  or  discontinuous  surface  slope 
onto  the  next  chordwise  region.  These  values  are  therefore 
used  only  on  regions  ending  in  the  interior  of  a  section.  The 
last  point  on  a  section  is  specified  by  NODEC  =  3  and  is  the 
only  node  point  that  must  always  be  specified  even  if  manual 
panelling  has  been  selected.  Negative  NODEC  values  are  also  per¬ 
mitted  and  initiate  a  special  copying  routine  described  in  3.2.5. 

Four  panel  spacing  options  are  provided  in  the  A.P.R.  The 
action  of  INTC  values  of  0,  1  and  2  is  illustrated  in  Figure  10 
and  is  based  on  the  cosine  distribution  giving  increased  panel 
density  towards,  respectively,  the  beginning  and  end,  the  begin¬ 
ning  only,  or  the  end  only,  of  the  region.  Equal  spacing  through¬ 
out  the  region  is  provided  by  ISPAC  =  3.  Coupled  with  the  flexi¬ 
bility  offered  by  the  choice  of  chordwise  region  location,  these 
spacing  options  have  proven  adequate  so  far?  however,  other 
options  could  easily  be  added  should  the  need  arise  later,  e.g., 
one  based  on  increments  in  integrated  surface  curvature,  or  on 
increments  in  doublet  value  from  a  preliminary  two-dimensional 
solution  for  the  section. 

Clearly,  node  cards  provide  the  user  with  an  extremely  versa¬ 
tile  panelling  tool.  With  one  card  deck  of  basic  points  defining 
the  configuration  geometry,  he  can,  from  run  to  run,  change  the 
form  of  the  panelling  simply  by  changing  two  integer  values  on 
each  node  card.  Not  only  that,  he  can  also  move  node  cards  within 
the  deck  (but  not  the  node  cards  at  section  ends)  or  remove  some 
or  add  new  ones  from  run  to  run.  This  allows  the  user  to  concen¬ 
trate  his  panelling  in  the  areas  of  interest,  leaving  other  areas 
more  sparsely  panelled.  It  thereby  provides  a  very  effective  use 
of  the  limited  number  of  panels  available,  yet,  on  a  subsequent 
run  a  few  small  changes  to  the  node  cards  allow  the  emphasis  to  be 
switched  to  another  area  without  having  to  punch  a  new  basic  geo¬ 
metry  card  deck. 
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*,  -  SO  -cos(i,))/2 

WHERE  8,  -  (I  -  1  )w/N  AND  N  IS  THE  NUM8ER  OF  INTERVALS  REQUIRED 

(a)  INTC  •  0 


s,  =  S(1  -COS(8,)) 
WHERE  8,  *  (I  -  1  )w/2N 

(b)  INTC  =  I 


*J  »  S  SIN(6|) 

WHERE  0, -(l-1)w/2N 

(«)  INTC  =  2 


Figure  10. 


Spacing  Options  0,  1  and  2  in  the  A.P.R. 


In  most  of  our  applications,  the  section  basic  points  are 
assembled  in  two  chordwise  regions.  The  regions  correspond  to  the 
section's  lower  and  upper  surfaces,  respectively  (Figure  7).  The 
leading-edge  point  separates  the  two  regions  and  is  identified 
by  a  node  card  just  after  it  with  NODEC  =  1  and  JNTC  =  2.  The 
last  basic  point  on  a  section  is  followed  by  a  node  card  with 
NODEC  =  3  and  INTC  =  1. 

There  is  just  one  important  ground  rule  for  the  use  of  node 
cards:  the  total  number  of  panels  (automatic?  and/or  manual)  on 

each  section  of  a  patch  must  be  the  same.  The  total  is,  in  fact, 
the  number  of  panel  rows,  NROW,  for  that  patch.  The  program 
monitors  the  number  of  panels  on  each  section  and  the  calcula¬ 
tions  are  terminated  with  an  error  message  should  the  user  make 
a  mistake.  Provided  this  ground  rule  is  satisfied,  it  is  not 
necessary  for  the  panel  distribution  to  be  the  same  from  section 
to  section--in  other  words,  the  number  of  chordwise  regions  and 
their  node  information  can  vary  from  section  to  section. 


3.2.4  Spanwise  Regions 

Sections  defined  within  each  patch  may  be  assembled  in  a 
number  of  SPANWISE  REGIONS  for  the  purpose  of  controlling  panel 
density  and  spacing  in  the  spanwise  direction.  In  forming  span- 
wise  regions,  sections  defined  by  the  user  take  on  a  similar 
role  to  that  of  basic  points  in  the  chordwise  regions.  Although 
the  options  available  for  the  spanwise  regions  are  essentially 
the  same  as  described  for  the  chordwise  regions  in  3.2.3,  the 
two  are  applied  completely  independently;  for  example,  the  user 
may  request  automatic  panelling  in  the  chordwise  direction  and 
manual  in  the  spanwise  direction.  As  in  the  case  of  chordwise 
regions,  spanwise  regions  are  used  only  as  an  input  convenience 
and  are  discarded  once  the  panelling  is  complete. 

Spanwise  regions  must  end  at  user-defined  sections,  called 
NODE  SECTIONS ( Figure  11).  These  usually  coincide  with  kinks  in 
the  spanwise  direction  on  the  patch  planform,  but  the  user  can 
place  one  whenever  he  wishes  to  change  the  form  of  the  panelling 
or  to  change  between  manual  and  automatic  panelling  in  the  span- 
wise  direction.  For  convenience,  the  spanwise  node  information 
is  included  on  the  section  card  together  with  the  section  trans¬ 
formation  information  (3.2.2).  The  function  of  the  spanwise 
region  node  quantities,  NODS,  NPS ,  INTS--distinguished  from  the 
corresponding  chordwise  quantities  by  ending  in  S--follows  closely 
the  description  in  3.2.3.  NODS,  however,  must  be  set  to  zero 
(blank)  on  the  first  section  of  a  patch  and  on  all  intermediate 
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input  sections  that  are  not  node  sections.  (NPS  and  INTS  are 
then  inactive.)  The  last  section  on  a  patch  is  identified  by 
a  NODS  value  of  3  or  5;  5  is  used  if  the  patch  is  the  last  one 
on  the  configuration,  in  which  case  the  present  section  completes 
the  basic  description  of  surface  geometry. 

The  total  number  of  panels  defined  (manually  or  automatically) 
across  each  patch  in  the  spanwise  direction  is  monitored  by  the 
proqram  and  becomes  the  number  of  panel  columns,  NCOL,  for  that 
patch.  In  view  of  the  ease  of  qeneratinq  panels,  the  code  also 
monitors  the  runninq  total  of  panels,  and  if  a  limit  is  exceeded, 
the  calculation  terminates  with  an  appropriate  error  message. 

The  limit  is  set  internally  by  the  storage  capacity,  but  the  user 
is  given  the  opportunity  to  override  that  value  with  his  own 
estimate  (NPMAX  on  CARD  3)  of  the  total  he  intends  to  use  for  that 
case.  In  the  event  of  an  input  error,  this  will  avoid  the  in¬ 
advertent  and  expensive  use  of,  say,  1,000  panels  when  the  user 
intended  using  only  100. 

In  the  earlier  hover  work,®  it  was  noted  that  errors  in  the 
surface  doublet  distribution  can  develop  if  a  vortex  passes  within 
one  panel  width  of  the  blade  surface.  With  the  present  panelling 
options  in  the  spanwise  direction,  problems  of  this  nature  can  be 
minimized  by  requesting  a  patch  break,  edge,  near  the  first  vortex 
passage  and  selecting  appropriate  panel  numbers  in  adjacent  patches. 
For  this  purpose,  it  is  sometimes  convenient  to  know  the  panel 
number  needed  for  a  given  spacing  option  to  obtain  a  required 
minimum  panel  width,  say  AS^j^.  The  required  number  of  panels, 
NPSreq,  for  the  spacing  options  are  as  follows. 


INTS  =  0: 


NPS 


REQ 


2  A  S  M  T  ^T 
-1  ,  -j  MIN 

cos  (1  -  — = - 
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INTS  =  3: 
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( 6d) 


3.3  SPECIAL  ROUTINES 

The  geometry  routines  described  above  may  be  applied  for 
the  complete  configuration;  however,  special  routines  have  been 
provided  to  reduce  user  input  and,  in  particular,  to  avoid  du¬ 
plicating  information  already  supplied.  These  routines,  which 
are  described  below,  are  optional. 


3.3.1  Copying  Routine 

We  have  already  seen  (3.2.2)  a  copying  facility  accessible 
at  the  section  input  level.  This  copies  over  a  complete  sec¬ 
tion,  including  the  chordwise  region  information,  and  has, 
therefore,  a  rather  limited  application.  More  general  copying 
routines  are  provided  and  are  activated  at  the  basic  point 
level  to  copy  STRINGS  OF  BASIC  POINTS,  rather  than  complete 
sections.  This  capability  allows  a  new  section  to  be  assembled 
from  parts  of  previously  defined  sections.  Several  strings  of 
basic  points  may  be  assembled  from  a  number  of  previously  defined 
sections  and  the  points  selected  need  not  follow  the  same  direc¬ 
tion  as  originally  specified.  Furthermore,  the  copied  strings  of 
points  may  be  intermixed  with  strings  of  manually  input  basic 
points  to  complete  the  new  section. 
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For  this  copying  mode,  the  value  of  INMODE  on  the  section 
card  (3.2.2)  must  be  in  the  range  1  to  4.  'The  copying  is  acti¬ 
vated  by  inserting  a  NODE  CARD  having  a  NEGATIVE  sign  on  NODEC . 
This  is  regarded  as  a  DUMMY  node  card  because  it  does  not  neces¬ 
sarily  terminate  a  chordwise  region  (see  below) .  The  negative 
value  for  NODEC  determines  the  action  at  the  end  of  the  copied 
string  of  basic  points.  If  NODEC  =  -1  or  -2,  then  the  last 
copied  point  becomes  the  end  of  a  chordwise  region  on  the  new 
section  and  signifies,  respectively,  continuous  or  discontinuous 
slope  onto  the  next  chordwise  region.  We  then  continue  to 
specify  further  basic  points  or,  by  inputting  another  negative 
node  card,  we  can  copy  another  string  of  basic  points,  and  so  on. 
If  NODEC  =  -3,  then  the  last  copied  point  in  the  string  completes 
the  new  section. 

If  the  user  does  not  require  a  chordwise  region  to  end  at 
the  last  point  in  a  copied  string,  then  he  sets  NODEC  =  -4 
when  he  initiates  the  copy.  When  the  string  has  been  copied 
over,  the  program  then  expects  to  receive  further  basic  points 
to  complete  the  chordwise  region  or  another  negative  node  card 
can  be  used  to  copy  another  string  of  points,  and  so  on.  Clearly, 
if  NODEC  =  -4,  then  the  NPC  and  INTC  values  on  the  NODE  CARD 
are  inactive  and  may  be  left  blank. 

Whenever  a  negative  node  card  is  inserted,  it  must  always 
be  followed  by  a  COPY  CARD  containing  the  following  information 
(four  integers)  defining  the  location  of  the  required  string 
of  points,  IPCH,  ISEC ,  IB,  LB. 

IPCH  is  the  patch  number  containing  the  reauired  points. 

ISEC  is  the  section  number  relative  to  the  start  of 

that  patch. 

IB,  LB  are,  respectively,  the  first  and  last  basic  point 
numbers  (inclusive)  defininq  the  string.  The 
numbering  is  relative  to  the  start  of  the  section 
ISEC. 

Thus,  even  in  a  complicated  configuration,  it  is  relatively 
easy  to  specify  a  string  of  basic  points. 

This  option  offers  not  only  an  alternative  to  the  earlier 
copying  routine,  but  also  a  more  general  capability  because  the 
copying  is  initiated  at  the  basic  point  input  level,  rather  than 
at  the  section  input  level.  For  example,  the  complete  copied 
section  need  form  only  a  part  of  the  new  section,  it  being 
possible  to  have  other  basic  points,  both  before  and  after  the 
copied  string.  In  addition  to  this,  the  ability  to  break  the 
copying  into  strings  of  points  allows  a  new  distribution  of  chord- 
wise  regions  to  be  selected. 
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One  restriction  must  be  considered  when  using  this  copy 
routine--the  new  secticv  s  value  for  INMODE  must  coincide  with 
the  INMODE  values  on  s  ction  from  which  strings  of  points  are 
to  be  copied.  This  restriction  has  not  posed  a  problem  so  far, 
but  if  it  does,  it  would  not  be  too  difficult  to  remove. 


3.3.2  Auto  atic  Patch  Generator 

The  patch  covering  the  tip  edge  can  be  input  by  the  user 
as  an  ordinary  'itch,  but  this  can  get  tedious.  Optional  auto¬ 
matic  procedure  have  bee.*  installed  which  simplify  this  input 
bv  generatina  ^omelet e  patch  within  the  code.  This  AUTOMATIC 

PATCH  GENERATOR  A.P.G.)  is  initiated  at  the  patch  incut  level  bv 
inserting  a  nonzei.  vaiao  ror  parameter,  MAKE,  on  the  patch 
data  card.  The  value  '-•f  MAKE  identifies  the  patch  number  on 
the  edqe  cf  which  a  closing  patch  is  to  be  generated.  The  sign  of 
MAKE  determines  whether  thv  new  patch  is  on  side  3  (positive)  or 
side  1  (negative)  of  the  basic  patch. 

Consider  Figure  12.  Here  we  have  already  defined  the  patch 
representing  the  main  surface.  The  end  section  of  that  patch 
provides  the  BASE  SECTION  from  which  the  A.P.G.  creates  the  new 
patch  according  to  user  instructions.  When  the  A.P.G.  has  been 
activated,  the  next  card  must  contain  the  following: 

NPC,  INTC  ,  KURV ,  NODS,  NPS  ,  INTS . 

The  generated  patch  has  one  "chordwise”  region  with  NPC  panels 
spaced  according  to  the  value  of  INTC.  The  values  for  NPS  and 
INTS  determine  the  "spanwise"  panelling  on  the  patch.  It  is 
recommended  that  the  same  values  be  used  for  these  as  were  used, 
respectively,  for  NPC,  INTC  on  the  lower  surface  of  the  base  sec¬ 
tion;  In  particular,  NPS  =  0  should  only  be  used  if  the  value 
for  NPC  on  the  base  section  was  also  zero.  The  value  of  NODS  must 
be  5  as  this  patch  completes  the  patch  input. 

The  function  of  KURV  is  described  next.  Sections  defining 
the  new  patch  are  created  automatically  from  the  base  section 
coordinates.  The  contour  of  each  section  generated  may  be  either 
a  straight  line  ("square-cut”  tip)  or  a  semicircle,  depending  on 
the  value  of  the  quantity,  KURV,  supplied  by  the  user.  If  KURV 
is  0,  sets  of  basic  points  are  generated  on  straiqht  lines  join- 
inq  upper  and  lower  points  on  the  base  section.  The  same  number 
of  points  is  created  even  if  the  interval  across  the  base  section 
is  zero ^  e.g.,  at  the  leading  and  trailing  edges  (Figure  12). 

If  KURV  is  1,  the  basic  points  are  created  on  semicircles 
having  a  diameter  equal  to  the  local  "thickness”  of  the  base 
section . 


(•)  PUT  EDOI  PATCH  (KURV  -  0) 


(b)  PATCH  WITH  SEMICIRCULAR  SECTIONS  (KURV  =  1) 


Fiqure  12.  Automatic  Patch  Generator. 
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4.0  WAKE  ROUTINES 


In  general,  the  prescribed  wake  methods  and  the  relaxed  wake 
techniques  developed  for  the  HOVER  code  have  been  carried  over 
for  ROTAIR.  Basic  changes  that  were  required  for  the  surface 
potential  method  involved  the  attachment  of  the  wake  to  the 
actual  blade  surface  at  its  coned  and  twisted  position,  the  struc¬ 
tural  change  for  the  wake  from  individual  filaments  to  wake 
columns  for  the  velocity  potential  calculation,  the  description 
of  the  rotor  wake  in  patches  of  panels,  the  modelling  of  the 
initial  tip-edge  wake  geometry  that  separates  from  the  thick-blade 
surface,  and  a  velocity  potential  model  of  the  far-wake  influences. 


4.1  WAKE  PATCHING 


Within  ROTAIR,  the  rotor  wake  structure  is  now  actually 
composed  of  three  wakes  as  illustrated  in  Figure  13.  WAKE  1  con¬ 
sists  of  all  the  wake  columns  that  make  up  the  inner  sheet,  i.e., 
the  shed  wake  up  to  the  maximum  circulation  location.  WAKE  2 
is  composed  of  the  wake  columns  comprising  the  outer  tip  sheet  to 
the  trailing  edge  of  the  tip  section;  therefore,  the  first  wake 
column  of  WAKE  2  is  the  potential  representation  of  the  merged 
tip  vortex.  WAKE  3  includes  the  columns  of  wake  panels  that  are 
shed  from  the  tip-edge  separation  line.  The  extent  of  the  sepa¬ 
ration  across  the  blade  tip  is  controlled  by  the  user  and  is 
described  later  in  this  section.  Regardless,  WAKE  1  and  WAKE  2 
are  always  present  while  WAKE  3  can  be  eliminated  if  no  tip-edge 
separation  is  required. 

As  shown  schematically  in  Figure  13,  the  attachment  of  the 
wake  patches  to  blade  edges  is  independent  of  the  patches  of  panels 
representing  the  blade  surface.  For  example,  the  attachment  line 
of  WAKE  1  can  terminate  in  the  middle  of  a  blade  surface  patch. 
Also,  the  actual  streamwise  lengths  of  each  wake  are  controlled 
independently  in  the  manner  set  up  in  the  HOVER  code.  Of  course, 
the  length  of  WAKE  3,  if  it  exists,  is  set  by  the  user-specified 
merger  azimuthal  angle,  PSIMRG. 

In  order  to  simplify  user  input,  the  actual  patching  of  the 
wakes  and  the  stitching  of  the  wakes  to  the  blade  surface  are 
automatically  computed  in  ROTAIR.  In  this  way,  user  input  of  a 
rotor  wake  for  ROTAIR  is  almost  identical  to  that  required  for 
HOVER. 

4.2  INITIAL  WAKE  GEOMETRY 


The  global  wake  structure  is  still  represented  by  the  near-, 
intermediate-  and  far-wake  regions  originally  developed  for  the 
HOVER  code.  However,  the  new  far-wake  model  described  later 
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13.  Wake  Patching. 
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is  now  keyed  to  a  PSIF2  specification  so  that  the  far-wake 
boundary  is  closer  to  the  blade  for  inner  sheet  filaments  (Figure 
14) -  This  change  has  resulted  in  significant  savings  of  computer 
time;  however,  the  old  model  may  be  required  for  wake  relaxation 
calculations. 

The  near  wake  normally  includes  a  few  vortex  passes  below 
a  blade,  and  the  general  wake  equations  (3)  describing  the  fila¬ 
ments  in  this  region  are  summarized  below  (see  Figure  15). 


•  2  TT 

(1)  Tip  Vortex  Axial  Coordinates:  =  -g- 
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(2)  Tip  Vortex  Radial  Coordinates: 
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(3)  Vortex  Sheet  Axial  Coordinates;  (usually,  ipQ  =  j) 
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(4)  Vortex  Sheet  Radial  Coordinates: 
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As  shown  in  Figure  15,  zt  are  the  axial  coordinates  with  respect 
to  the  tip  path  plane.  However,  the  wake  actually  leaves  from 
the  blade  trailing-edge  position.  Consequently,  an  effective  ki 
is  calculated  in  ROTAIR  to  carry  the  tip  wake  filament  from  the 
trailing  edge  to  the  axial  position  at  first  blade  passage.  The 
equation  for  the  tip  vortex  with  respect  to  the  tip  path  plane 
system  is  then 
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equation  is  similarly  modified. 


These  wake  coordinates  must  be  transformed  to  the  H.C.S. 
for  use  with  the  blade  surface  panelling.  In  general,  the  blade 
cone  angle,  3,  requires  a  vertical  displacement  and  a  small  in¬ 
board  shift  of  the  tip  coordinates  as  shown  in  Figure  16. 


ZHCS'  2 


LHCS 


Figure  16.  Effect  of  Blade  Coning  on  Wake  Coordinates. 


Now,  Az  =  (R  -  r^)  +  sin  3 

and  Ar  =  R  -  £ (R  -  rf)  cos  3  + 

Hence , 

ra  =  R  -  Ar  =  (R  -  rf)  cos  3  +  rf  . 

Consequently,  the  tip  path  plane  tip  vortex  coordinates  (Eq 
(7a),  (7b)  and  (7c))  are  modified  to  obtain  the  H.C.S.  coordinates 

as  follows: 


HCS  =  Zt  +  Az 

(8a) 

HCS  r  ra 

(8b) 
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The  detailed  near  wake  segmentation  is  modelled  in  the  manner 
of  the  HOVER  code  (Figures  17  and  18).  The  azimuthal  step  angle, 
Aij^,  of  the  wake  sheet  and  the  tip  region  can  be  specified  indepen¬ 
dently,  i.e..  Ail  and  A\p  .  Further,  these  step  sizes  in  the 

s  wt 

wake  (or,  equivalently,  the  vortex  segment  lengths)  can  be  changed 
to  new  values  after  the  first  blade  passage.  This  flexibility 
allows  very  efficient  computation  without  sacrificing  accuracy 
in  loads  prediction. 

The  tip  separation  wake,  WAKE  3,  is  modelled  as  illustrated 
in  Figure  19.  The  user  specifies  the  number  of  surface  panels 
from  the  trailing  edge  that  participate  in  the  wake  separation, 
NSEP;  the  height  above  the  trailing  edge  that  the  tip  vortex 
attains,  DELZ;  and  the  inboard  shift  of  the  tip  vortex  across  the 
blade  surface  at  the  trailing  edge,  DELR.  As  a  first  approxima¬ 
tion,  DELZ  and  DELR  can  be  set  to  a  distance  corresponding  to  the 
tip  vortex  leaving  the  blade  surface  at  1/2  the  tip  section  geo¬ 
metric  angle  of  attack.  As  experience  is  gained  with  the  tip 
shedding,  this  procedure  will  be  automated  within  the  code. 

For  this  model,  the  feeding  sheet  filament  effects  are  in¬ 
cluded,  and  the  actual  wake  panelling  of  WAKE  3  is  automatically 
generated. 

The  intermediate-wake  region  is  composed  of  fixed-pitch, 
constant  radius  helices.  The  pitch  and  radius  of  each  filament 
is  determined  by  the  final  pitch  and  radius  of  the  filament  in 
the  near  wake. 

Finally,  recall  that  in  the  HOVER  code  the  far-wake  region 
represents  a  semi-infinite  continuation  of  this  vorticity  model. 
That  is,  the  downwash  velocity  components  due  to  each  helical 
vortex  filament  trailing  behind  the  rotor  is  approximated  by  the 
velocity  due  to  a  cylinder  of  uniform  vorticity  in  the  far  wake 
(Figure  20) . 
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Figure  20,  Far-Wake  Vorticity  Model. 


Further,  the  velocity  due  to  a  semi-infinite  cylinder  of 
uniform  vorticity  is  identical  to  that  induced  by  a  sink  sheet 
on  the  capping  surface  of  uniform  strength  equal  to  the  vorticity 
( figure  21). 

In  ROTAIR,  this  same  model  is  used  for  velocity  calculations 
with  the  modification  that  the  starting  axial  location  of  the  far- 
wake  influence  is  individually  determined  for  each  vortex  helix. 
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piqure  21.  Equivalency  of  a  Cylinder  of  Uniform  Vorticitv  and 
a  Sink  Sheet  of  Uniform  Strength. 

While  this  sink  sheet  method  works  well  for  velocities,  such 
a  simplification  for  the  velocity  potential  influence  is  not  pos¬ 
sible.  Instead,  the  velocity  potential  induced  by  the  far-wake 
is  that  due  to  the  superposition  of  cylindrical  potential  sheaths 
of  linearly  increasing  doublet  distribution.  The  velocity  and 
potential  equivalence  is  schematically  shown  in  the  figure  below. 


Velocity  Potential 


Figure  22.  Illustration  of  the  Velocity  Potential 

Produced  by  a  Cylinder  of  Uniform  \orticity. 
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This  far-wake  potential  model  has  been  added  to  ROTAIR  for  the 
aerodynamic  matrix  calculation.  The  accuracy  of  the  model  is 
demonstrated  in  Section  6.0. 


4 . 3  PRESCRIBED  WAKE  ITERATION 

The  prescribed  wake  constants  are  determined  by  the  empirical 
equations  due  to  Landgrebe  3  or  Kocurek.  4  These  equations 
are  also  described  in  detail  in  the  HOVER  document  .  ^  Regard¬ 
less  of  which  equations  are  used,  the  wake  constants  are  related 
to  thrust  coefficient  (more  recently,  the  constants  have  been 
related  to  maximum  blade  circulation,  and  not  to  blade  collective). 
Therefore,  the  prescribed  wake  iteration  should  be  orooerly 
carried  out  within  a  collective  iteration.  Consequently,  for  this 
intermediate  time  until  the  collective  loop  is  available  in  ROTAIR, 
the  user  is  required  to  read  in  the  appropriate  wake  constants 
explicitly.  Several  computer  runs  will  then  be  required  to  obtain 
the  collective  setting  for  a  calculated  thrust  coefficient  equal 
to  the  requested  thrust  coefficient  used  to  position  the  wake. 

The  prescribed  wake  iteration  loop  in  ROTAIR  continues  until  the 
blade  maximum  circulation  location  is  unchanged. 


4. 4  RELAXED  WAKE 

If  requested,  the  merged  tip  vortex  position  is  then  relaxed 
to  a  force-free  location  for  the  resultant  fixed  collective.  The 
relaxation  technologies  developed  for  HOVER  are  incorporated  here. 
These  were  all  described  in  great  detail  in  Reference  8,  and  the 
reader  is  referred  to  that  document. 


kocurek,  J.D.,and  Tangier,  J.L.,  "Prescribed  Wake  Lifting 
Surface  Hover  Performance  Analysis",  Presented  at  the  32nd 
Annual  National  VSTOL  Forum  of  the  American  Helicopter 
Society,  Washington,  D.C. ,  May  1976. 
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5.0  ROTOR  PERFORMANCE  CALCULATIONS 

Since  ROTAIR  is  a  surface  singularity  method,  the  local 
surface  velocities  are  calculated  at  all  blade  surface  panel 
center  points.  With  this  information,  the  local  panel  surface 
pressure  coefficients  are  obtained  and  rotor  loading  can  then  be 
calculated  by  surface  pressure  integration.  Usually,  the  pres¬ 
sure  coefficient  is  defined  by  reference  to  the  local  sectional 
onset  flow  velocity,  say  VQ.  For  the  rotor  problem,  this  quantity 
varies  linearly  with  the  blade  radius  station.  For  the  hover 
problem  considered  here,  the  characteristic  velocity  for  normali¬ 
zation  is  the  rotor  tip  rotation  speed,  Vtip  -  ftR.  Consequently, 
it  is  convenient  to  define  two  pressure  coef f icients--a  conven¬ 
tional  one  based  on  the  local  sectional  onset  velocity  for  com¬ 
parison  with  experimental  data,  and  another  based  on  the  tip 
rotation  velocity  to  observe  the  actual  blade  surface  pressure 
distribution  and  also  to  facilitate  the  direct  calculation  of 
the  integrated  rotor  performance.  As  indicated  in  Figure  23, 
the  sectional  characteristic  onset  flow  velocity  at  section  i 
is  Vc^ ,  so  that 


V 

o  . 
1 


WCLIMB 


+  ft2Y  .  2 
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Consequently,  the  conventional  pressure  coefficient  at  panel  j 
of  section  i  in  variables  normalized  by  rotor  radius,  R,  is 


C 


w  2  +  y  .  2 
c  y  i 


(9) 


On  the  other  hand,  the  pressure  coefficient  at  panel  i  based  on 
rotation  tip  speed  is 
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(10) 
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figure  23.  Sectional  Aerodynamic  forces. 
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In  both  equations,  r\^  is  the  normalized  radius  to  panel  j,  and 

v  is  the  normalized  surface  velocity  at  panel  j. 
s  .  . 

ID 


5. 1  INVISCID  ROTOR  AERODYNAMIC  PROPERTIES 

The  sectional  properties  can  be  found  by  integration  of  the 
pressure  distribution  around  the  airfoil  surface.  For  our  pur¬ 
poses  here,  the  sectional  loads  are  calculated  for  a  given  column 
of  panels.  The  orientation  of  the  column  may  or  may  not  be 
streamwise,  depending  on  the  user  input  surface  panelling  requests. 
The  following  figure  will  be  useful  in  the  derivation  of  the 
sectional  loading  (see  also  Figure  23) . 


Figure  24.  Illustration  of  the  Pressure  Loading  on 
Panel  j  of  Column  i. 


By  definition 
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Integrated  Induced  Torque  Coefficient: 


•  k 


The  sectional  force  and  moment  components  along  the  other  H.C.S. 
axis  are  also  computed  in  the  program.  Further,  the  sectional 
lift  and  induced  drag  coefficients  are  obtained  from  Eq  (11)  in 
the  manner  described  in  References.  That  is  (see  Figure  23), 
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Finally,  if  the  mean  velocity  at  a  section  is  approximated 
by  the  section  onset  flow  only,  then  the  sectional  lift,  thrust, 
and  integrated  thrust  coefficients  can  be  obtained  from  the 
sectional  circulation  value,  r^,  as  was  done  in  Reference  8;  i.e., 
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These  exoressions  are  useful  in  that  normal  loads  calculated  by 
circulation  converge  faster,  i.e.,  with  fewer  surface  panels 
than  those  calculated  by  pressure  integration. 


5. 2  VISCOUS  AERODYNAMIC  CONTRIBUTIONS 

ROTAIR  includes  the  capability  of  calculating  the  profile 
losses  in  the  same  manner  as  the  HOVER  code.  However,  in  ROTAIR 
the  defining  stations  for  two-dimensional  aerodynamic  tables  are 
independent  of  the  stations  used  to  describe  the  blade  geometry. 
This  greatly  simplifies  the  user  input  of  the  aerodynamic  tables 
for  complex  rotors.  The  profile  drag  coefficients  are  obtained 
by  table  interpolation  at  the  calculated  section  lift  coeffi¬ 
cients  and  the  resultant  profile  torque  and  thrust  losses  are 
computed.  The  reader  is  referred  to  Reference  8  for  details. 
Eventually,  this  two-dimensional  empirical  calculation  of  profile 
losses  should  be  replaced  with  a  prooer  boundary  layer  calcula¬ 
tion  since  the  actual  surface  pressure  distributions  are  now 
available . 
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6 . 0  PRELIMINARY  RESULTS 


The  basic  doublet  code  (VSAERO)  for  fixed-wing  analysis^ 
has  been  checked  for  many  different  wing  configurations  includ¬ 
ing  the  detailed  flow  around  wing  tips.  Still,  this  code  is 
constantly  being  improved  and  validated.  In  addition,  the  rotor 
version  of  the  code,  ROTAIR,  described  in  this  report,  has  been 
only  recently  developed  and  has  been  used  for  only  a  few  con¬ 
ventional  rotor  configurations.  Consequently,  the  results  pre¬ 
sented  here  are  preliminary.  Example  calculations  presented 
below  serve  to  validate  the  programming  and  to  illustrate  the 
capabilities  and  limitations  of  ROTAIR. 


6. 1  FIXED-WING  TIP  FLOWS 

Fixed-wing  calculations  have  demonstrated  the  potential  of 
the  code  for  the  detailed  analysis  of  flows  around  rotor  tips.9 
For  example,  steady  pressure  distributions  calculated  on  a  rec¬ 
tangular  wing  and  plotted  spanwise  against  surface  distance  in 
a  number  of  cross-flow  planes  show  a  steady  growth  of  the  peak 
suction  associated  with  attached  flow  (Figure  25) .  The  wing 
has  a  rounded  tip  surface.  The  suction  peak,  which  dominates  the 
distributions  downstream  of  the  half-chord  location,  arises  from 
growing  velocity  components  in  the  cross-flow  planes.  Favorable 
pressure  gradients  exist  in  cross-flow  planes  near  the  leading 
edge,  but  adverse  gradients  develop  once  the  suction  peak  becomes 
established,  and  these  eventually  reach  unrealistic  values  under 
the  attached  flow  constraints. 

Prescribing  separation  along  the  tip  does  lead  to  more 
realistic  tip  pressure  distributions.  However,  comparisons  with 
experimental  measurements  have  indicated  the  need  to  establish 
the  extent  of  the  edge  separation  in  a  more  consistent  manner. 

For  example,  calculated  pressures  on  a  swept  wing  are  shown  in 
Figure  26  plotted  normal  to  the  surface  in  a  cut  through  the 
tip  at  .72  x/c.  Incidence  is  13°  and  the  wing  is  twisted  4° 
nose  down  at  the  tip  about  the  half-chord  line.  The  extent  of 
edge  separation  was  set  at  about  30%  of  the  tip  chord  in  this 
case.  The  calculations  for  three  wake  iteration  cycles  are 
compared  with  experimental  measurements  taken  by  the  N.L.R. 

The  experimental  values  have  been  interpolated  at  the  calcula¬ 
tion  station.  The  general  agreement  appears  good,  but  the 
upoer  surface  suction  is  underpredicted  by  the  calculation. 
Although  the  spanwise  velocity  component  on  the  upper  surface  has 
reversed  in  the  calculation,  its  level  is  clearly  too  low. 

Also  the  tip  lower  surface  pressure  is  low  compared  with  the 


^5Labrujere,  Th.E.,and  De  Vries,  0.,  "Evaluation  of  a  Poten¬ 
tial  Theoretical  Model  for  the  Wake  Behind  a  Wing  via  Comparison 
of  Measurements  and  Calculations",  NLR ,  TR-7 4036  U,  July  1974. 
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Figure  25.  Spanwise  Pressure  Distributions  Near  the  Tip 
of  a  Rectangular  Wing  in  Steady  Flow. 
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EXPERIMENT  (REF.  8) 

PRESENT  CALCULATIONS 
WAKE  ITERATION 


26.  Snanwise  Pressure  Distribution  Calculated  Near  the 
Tin  of  a  Swept  Winq  with  Prescribed  Edqe  Seoaration 


experimental  value,  indicating  that  the  cross-flow  velocities 
there  are  too  high  in  the  calculation.  Referr-^.g  to  the  cross- 
flow  analoqy,  it  would  appear  that  these  discrepancies  are  con¬ 
sistent  with  too  early  a  time;  i.e.,  the  calculation  needs  a 
larqer  extent  of  separation  in  this  case. 

Further  work  is  planned  at  AMI  to  investigate  ways  of 
establishing  the  extent  of  separation  along  the  tip  edge.  The 
couolinq  of  a  boundary  layer  calculation  will  be  considered. 
Meanwhile,  the  present  method  provides  a  tool  for  the  investiga¬ 
tion  of  the  tip  flow  and  vortex  formation. 


6.2  ROTOR  THICKNESS  EFFECTS 

A  sample  single-bladed  rotor  case  was  used  to  investigate 
the  effects  of  blade  thickness.  In  order  to  eliminate  the 
effects  of  a  returning  wake,  the  wake  was  truncated  at  180° 

(h  revolution)  and  the  far-wake  influence  was  set  to  zero.  The 
blade  was  an  aspect  ratio  16,  untwisted  rectangular  rotor. 
Compressibi li ty  was  eliminated  and  the  rotor  was  set  to  1°  of 
collective  in  order  to  maximize  the  thickness  percentage 
changes.  The  results  are  illustrated  in  Figure  27.  As  can  be 
seen  on  this  expanded  scale ,  the  2%  thick  airfoil  results  com¬ 
pare  very  well  with  the  zero  thickness  vortex- lattice  results. 
Increasing  the  thickness  to  12%  results  in  a  loss  of  circulation. 
This  trend  is  just  the  opposite  of  what  occurs  in  the  two-dimen¬ 
sional  case.  An  empirical  equation  that  is  commonly  used  for 


wing  design  is 

/dCl\ 

/  t  \ 

\  da /THICK 

(1  +  f5) 

\  da  /  tHIN  ' 

where  f  '  0.8.  A  two-dimensional  version  of  the  surface  singu¬ 

larity  code  and  VSAERO^  were  used  to  examine  the  thickness  effects 
for  fixed  wings.  The  following  results  were  obtained. 
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The  two-dimensional  results  were  for  a  Joukowski  airfoil.  The 
VS AERO  runs  were  for  t/c  =  0.12  and  t/c  =  0.02.  The  results  for 
the  aspect  ratio  six  wing  included  increased  tip  loading  to 
simulate  rotor  loading  by  linearly  twisting  the  wing  tip  upward 
8°.  Evidently,  thickness  effects  are  reversed,  or,  at  least 
greatly  reduced  in  the  presence  of  highly  three-dimenisonal 
flows.  This  behavior  should  be  more  thoroughly  investigated  in 
the  future. 

Similar  comparison  calculations  were  made  for  a  two-bladed 
rotor  at  1°  collective  with  fixed-wake  geometry.  These  are 
shown  in  Figure  28.  Once  more,  thickness  reduces  the  blade 
loading  for  the  two-bladed  case. 


6 . 3  A  CLOSE  VORTEX  PASSAGE  CASE 

Figures  29  through  31  illustrate  various  calculated  per¬ 
formance  characteristics  for  a  single-bladed ,  untwisted  rotor 
with  an  NACA  0012  airfoil  section.  As  this  was  a  preliminary 
calculation,  the  tip  vortex  shedding  was  not  modelled  (vortex 
wake  shedding  only  from  the  trailing  edge).  Also,  only  ten 
panels  were  used  around  the  airfoil  section  and  ten  columns  of 
panels  along  the  blade  radius.  While  this  is  too  few  panels  for 
obtaining  all  of  the  performance  details  in  the  neighborhood  of 
the  vortex,  it  is  enough  to  verify  the  calculation  capabilities 
of  the  code. 

Two  different  calculations  were  made  for  comparison.  The 
basic  tip  vortex  wake  constants  for  the  two  cases  were 
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The  radial  position  of  the  tip  vortex  for  both  cases  was 
=  0.875,  However,  the  first  blade  passage  axial  distance 
for  case  "a"  was  z^  -  -0.0377,  while  that  for  case  "b"  was 

z*b  ”  -0,126.  The  blade  trailing  edge  was  at  zTE  =“0.0125.  For 

these  calculations  the  wakes  were  arbitrarily  truncated  at  i{' 

=  480°.  w 

In  all  of  the  figures,  the  "dashed1’  lines  are  values  along 
the  blade  lower  surface  and  "solid"  lines  are  values  along  the 
blade  upper  surface.  Also,  the  blade  surface  geometries  at  the 
requested  data  cuts  are  overlaid  on  the  figures.  Figures  "a" 
illustrate  results  for  the  close  vortex  passage,  while 
"b"  show  results  for  the  distant  vortex  passaqe. 

The  radial  distribution  of  doublet  strength  is  shown  in 
Figure  29.  By  comparing  (a)  and  (b)  ,  the  buildup  and  the 
spreading  of  the  vorticity  near  the  vortex  and  across  the  blade 
chord  may  be  seen. 

The  strong,  direct  vortex  action  is  almost  wholly  confined 
to  the  blade  lower  surface,  while  the  upper  surface  values  are 
the  secondary  result  of  the  vortex  induced  doublet  distribution 
on  the  blade  surface.  This  represents  a  "classic"  vortex  surface 
interaction  effect  on  the  doublet  distribution  near  the  trailing 
edge  of  the  thick  blade  (recall  that  vorticity  is  the  gradient  of 
the  doublet  distribution) . 

Figure  30  compares  the  spanwise  velocity  component  varia¬ 
tion  near  the  blade  leading  edge.  The  strong  deceleration 
caused  by  the  free  vortex  to  the  usual  acceleration  of  the  flow 
on  the  lower  surface  around  the  tip  is  seen  by  comparing  the 
results  in  figures  "a"  and  "b".  Finally,  Figure  31  illustrates 
the  radial  distribution  of  minus  Cp  for  this  same  location  near 
the  blade  leading  edge. 

The  chordwise  pressure  plots  (not  shown  here)  display  the 
expected  reduction  in  the  suction  pressure  peak  inboard  of  the 
vortex  and  augmentation  of  the  suction  peak  outboard. 


6.4  OH- 5 8 A  ROTOR  C ALCULATIONS 


ratio 
per  fo 


The  OH-58A  helicopter  rotor  is  a  two-bladed  high  aspect 
rotor  with  -10.6  deqrees  of  linear  twist.  Extensive 
rmancc  prediction  comparisons  for  the  OH-58A  bv  the 
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figure  31.  Concluded 
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UTRC  and  HOVER  programs  have  already  been  reported  in  Reference  8. 
Calculations  by  these  programs  have  proven  very  useful  in  validating 
ROTAIR  in  examining  the  convergence  behavior  of  ROT AIR  with  panel 
number. 

The  blade  surface  panelling  of  this  rotor  is  shown  in  Figure 
32(a)  and  (b) .  In  this  case,  the  rotor  collective  angle  is  5.75°, 
and  the  blade  coning  angle  is  3°  to  correspond  to  the  previous 
calculations.  Additionally,  the  radial  panelling  distribution 
(15  columns)  is  equivalent  to  the  panelling  of  the  simpler  methods, 
and  the  Landqrebe^  wake  constants  were  set  for  a  thrust  coefficient 
of  0.0022.  Therefore,  the  wake  geometry  should  be  approximately 
eaual  for  the  three  methods. 


6.4.1  Far-Wake  Convergence 


ROTAIR  was  first  run  with  different  wake  lengths  to 
determine  the  convergence  behavior  of  the  new  far-wake  velocity 
potential  model.  With  only  ten  panels  around  the  airfoil  (i.e., 
10  rows)  and  30°  wake  segments,  the  convergence  of  several  key 
performance  properties  was  as  follows: 
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These  changes  are  mostly  due  to  the  loss  in  ultimate  wake  con¬ 
traction  as  is  shortened.  Apparently,  the  far-wake  model 
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could  almost  be  taken  right  up  to  the  blade  surface  for 
prescribed  wake  calculations!  This  also  indicates  that  for  most 
of  the  wake,  a  relaxed  wake  scheme  based  on  sheaths  of  vorticity 
(iterating  on  vorticity)  miqht  be  possible.  At  any  rate,  for 
further  calculations,  the  wake  length  could  be  set  for  if/,.  =  480° 

and  vjjf  =  520°  with  little  loss  in  accuracy. 
r2 


6.4.2  Circulation  Distribution  Comparisons 

The  circulation  distributions  computed  by  three  codes  are 
compared  in  Figure  33.  For  this  figure,  all  of  the  calculated 
circulation  data  are  scaled  by  the  maximum  circulation  predicted 
by  ROTAIR  (CIRC (MAX)  =  0.008672).  All  the  methods  agree  reason¬ 
ably  over  the  inner  portion  of  the  blade  (the  thick  blade  solu¬ 
tion  indicates  lower  circulation  values  aqain  (see  Section  6.2)). 
Finally,  as  reported  in  Reference  8,  the  lifting  line  program, 
UTRC ,  overpredicts  the  circulation  values  in  the  tip  region 
compared  to  the  ROTAIR  and  HOVER  results  which  agree  favorably. 


6.4.3  Convergence  with  Panel  Number 

The  convergence  behavior  of  integrated  performance  parameters 
obtained  by  integration  of  the  surface  pressures  is  much  too  slow 
in  the  current  code.  Figure  34  illustrates  the  variation  of  the 
thrust  coefficient  computed  by  surface  pressures  and  computed  by 
circulation  values  for  chordwise  panel  numbers  (2  x  NPC)  ranging 
From  10  to  30.  In  all  these  calculations,  =  480°  and  il’  =  520°. 

rl  f2 

As  usual,  the  circulation  computed  value  converges  from  below, 
while  the  pressure  computed  value  converges  from  above.  Also,  the 
circulation  based  thrust  converges  much  faster,  but  is  more  re¬ 
stricted  in  its  use.  That  is,  regions  of  high  radial  velocities 
or  arbitrary  geometries  invalidate  the  circulation  base  thrust. 

In  fact,  even  for  this  case  with  Z,  0.03,  the  two  thrust  cal- 

culations  seem  to  be  converging  on  slightly  different  asymptotic 
values . 
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ROTAIR 


Calculated  Circulation  Distribution  Comparisons  for  the  OH-58A  Rotor 
(Normalized  by  the  Maximum  Circulation  Value  Computed  by  ROTAIR) . 


NUMBER  OP  CHORDWISE  PANET.S  (2  x  N^C) 


piciuro  34, 


E  f f  ect 
Thrust 


of  Number  o\ 
Coe  f  f i cicn t . 


Chordwise  Panels  (NROW)  on 
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Figure  35  illustrates  the  convergence  behavior  of  several 
key  rotor  performance  quantities  for  chordwise  panel  number. 

In  this  figure,  each  of  the  ordinate  values  are  scaled  by  their 
respective  ultimate  values  achieved  with  30  chordwise  panels; 
hence,  all  the  quantities  approach  a  value  of  1.0  at  30  panels. 
The  values  for  these  quantities  with  30  panels  are 


CIRC (MAX) 
CT 

CTGAM 

CQ 

FM 


0.08853 

0.002242 

0.002160 

0.00009273 

0.554 


Obviously,  the  induced  torque,  or  section  drag  variation  with 
panel  number  is  quite  large,  and,  even  with  30  panels,  the 
error  in  the  induced  torque  is  not  acceptable  for  efficient 
computation. 


This  is  an  area  of  ROTAIR  that  needs  additional  work. 
Presently,  loadings  are  obtained  by  summing  C*  x  area  for  all 
of  the  panels.  This  can  be  improved  by  usingche  surface  arc 
length  around  each  given  section  along  with  a  higher-order 
interpolation  through  Cp  values. 

An  additional  method  to  obtain  a  highly  accelerated  con¬ 
vergence  of  the  sectional  properties  would  include  a  two- 
dimensional  program  in  ROTAIR.  Once  the  section  panelling  is 
complete,  the  two-dimensional  program  could  be  used  to  generate 
appropriate  factors  that  could  be  applied  later  for  the  three- 
dimensional  calculations.  For  example,  the  inviscid  two- 
dimensional  drag  should  be  zero  for  all  C^'s.  For  a  given  air¬ 
foil  panelling,  the  two-dimensional  results  would  provide  the 
ACD  that  should  be  subtracted,  section  by  section,  for  the 
three-dimensional  calculation.  The  more  panels  used,  the 
smaller  the  ACD's.  The  lift  values  could  be  corrected  in  a 
similar  manner.  Finally,  the  additional  computing  time  for  the 
two-dimensional  calculations  would  be  very  small. 


n 

B 


82 


6.4.4  Pressure  Distributions 


Although  the  completed  integrated  loads  by  pressure 
integration  converges  slowly  in  the  current  program,  the  actual 
surface  pressure  values  converge  rapidly,  even  for  as  few  as  a 
total  of  eight  panels  around  the  airfoil  (see  Ref.  11) ,  and 
detailed  pressure  distribution  calculations  will  prove  very  use¬ 
ful  in  the  analysis  of  the  newer  rotor  tip  shapes. 

Figures  36  (a) ,  (b) ,  (c)  and  (d)  demonstrate  the  calculated 

chordwise  pressure  distributions  for  the  30  panel  case  of 
Section  6.4.3  at  several  radial  locations.  The  airfoil  panelling 
shapes  (seen  in  the  H.c.S.)  are  superimposed  on  these  figures 
and  illustrate  the  linear  twist  of  the  blade  geometry  and  the 
blade  coning. 


6.5  AMES  MODEL  ROTOR 


The  ultimate  validation  of  the  calculated  pressure  distri¬ 
bution  is  a  direct  comparison  with  measured  surface  pressure 
data.  Here,  we  present  a  comparison  of  the  calculated  surface 
pressure  distribution  on  a  low  aspect  ratio  blade  with  those 
measured  by  Caradonna  and  Tung .  6 


The  actual  blade  surface,  including  blade-tip  closure,  is 
represented  by  quadrilateral  panels  as  shown  in  Figure  37.  To 
compare  with  the  detailed  measurements,  a  row  of  18  panels  was 
used  around  the  airfoils  and  15  columns  across  the  blade 
radius.  In  addition,  27  panels  were  employed  to  close  the  rotor 
tip  for  a  total  of  297  panels.  For  production  calculations, 
about  half  this  number,  or  150  panels,  are  sufficient.  The 
wake  constants  were  set  for  a  far-wake  radius,  R^AR,  equal  to 
0.88.  This  ultimate  contraction  radius  was  taken  from  the 
measured  wake  data. 

Figures  38  (a) ,  (b)  and  (c)  illustrate  the  surface  pressure 

comparisons  at  the  radial  stations  where  the  experimental  data 
was  collected  and  excellent  agreement  is  shown  at  all  stations. 
Furthermore,  the  cost  of  these  calculations  on  a  commerical 
CYBER  176  was  only  $25  (approximately  $12  on  the  Langley  175). 
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icjure  36.  Chordwisc  Pur fac^  "ressurc  Distributions  (,'linus  Cu)  for  the  OI1-58A 
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Figure  36.  Concluded 
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Figure  38.  Continued 
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Figure  38.  Concluded 


This  was  achieved  because  of  the  very  accurate  far-wake  doublet 
model.  In  fact,  since  most  all  of  the  contraction  occurs 
before  first  blade  passage  for  this  case,  only  270°  of  total 
detailed  wake  were  needed  for  this  two-bladed  rotor.  Figure 
39  illustrates  the  boundary  where  the  "far-wake"  analytical 
model  was  included. 


6.6  RELAXED  WAKE  CALCULATIONS 


A  relaxed  wake  calculation  in  ROTAIR  is  demonstrated  in 
Figure  40  for  the  OH-58A  blade.  Since  the  basic  free-wake 
methodology  has  been  carried  over  from  Reference  8,  the  con¬ 
vergence  of  the  axial  coordinates  should  follow  the  same 
approximate  trends  obtained  with  HOVER.  Of  course,  the  wake 
structure  is  not  exactly  eaual  for  the  two  codes  since  ROTAIR 
includes  exact  blade  twist  and  coning,  and  has  a  slightly  dif¬ 
ferent  far-wake  model  (see  Section  4).  The  calculations  in 
Figure  40  do  show  the  highly  convergent  behavior  found  with 
HOVER,  but  the  first  blade  passage  axial  distance,  Z,  ,  con¬ 


verges  at  a  much  slower  rate  than  the  HOVER  calculations. 
Detailed  velocity  scan  comparisons  for  the  two  programs  have 
demonstrated  that  ROTAIR  predicted  wake  induced  velocities  are 
correct  when  the  wake  vortex  strengths  are  equal.  For  this 
calculation,  NPC  =  5  so  that  only  10  panels  were  used  around 
the  airfoil.  As  shown  in  Section  6.4.3,  there  are  errors  in 
the  circulation  values  for  so  few  panels;  consequent ly ,  more 
chordwise  panels  should  always  be  used  if  relaxation  is  desired 
Additional  relaxation  runs  with  dense  panelling  should  be  under 
taken  to  complete  ly  val  idate  this  part  of  the  code.  Also,  the 
method  utilized  in  HOVER  for  storing  the  influence  coefficients 
for  the  riqid  wakes  needs  to  be  added  to  ROTAIR  to  reduce  re¬ 
laxation  costs. 
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7.0  SUMMARY  AND  FUTURE  DEVELOPMENT 

A  surface  singularity  potential  flow  code,  ROTAI R ,  has  been 
assembled  for  the  calculation  of  detailed  surface  pressures  on 
rotors  in  hover  or  climb.  This  method  is  basically  a  marriage 
of  two  codes  developed  at  AMI:  the  fixed-wing  surface  singular¬ 
ity  doublet  code^  and  the  rotor  lif tinq-sur f ace  code. 8  The 
program  includes  a  tip  vortex  separation  model.  Also,  the  rotor 
tip  surface  is  panelled  so  that  pressures  are  calculated  right 
around  the  tip-edge  surface. 

Preliminary  calculations  have  verified  the  capabilities 
of  the  program  for  computing  blade  surface  properties  in  the 
presence  of  a  close-vortex  passage.  Additionally,  calculated 
pressure  distributions  compare  favorably  with  experimental 
data  for  a  low  aspect  ratio  two-bladed  rotor,  and  the  cal¬ 
culated  circulation  distribution  is  consistent  with  that  com¬ 
puted  earlier  with  the  lifting-surface  code.  Finally,  the 
newly  developed  far-wake  doublet  model  promises  to  keep  computing 
ef forts  practical . 

Meanwhile,  further  development  of  the  methods  is  continuing. 
The  rotor  wake  relaxation  schemes  are  being  modified  to  allow 
for  free  contraction  boundaries  and  to  reduce  computing  effects. 
Elastic  modellinq  and  collective  iterations  are  also  being 
added  to  the  program. 

In  the  immediate  future,  thickness  effects  should  be 
investicrated  in  more  detail;  the  extent  of  tip-edge  separation 
and  tip  vortex  formation  needs  additional  work;  the  convergence 
of  integrated  pressure  data  with  panel  number  must  be  improved 
for  practical  computation;  influences  of  compressibility  on  the 
surface  sinaularity  method  should  be  examined;  and,  eventually, 
the  addition  of  a  boundary  layer  method  should  be  considered 
to  eliminate  the  two-dimensional  profile  approximations  and  to 
more  realistically  study  the  rotor  tip  problem. 
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APPENDIX  A 
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Figure  a-1.  Complete  Subroutine  Map 


Figure  A-2.  Prescribed  Wake  Loop. 
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Figure  A-3.  Relaxed  Wake  Loop. 


APPENDIX  B:  INPUT  DESCRIPTION 


The  input  is  in  five  main  parts:  Basic  Data,  Wake 
Blade  Orientation,  Geometry  and  Aerodynamic  Data  Tables. 
B— 1  shows  the  assembly  of  the  input  cards. 


BASIC  DATA 


CARD  1:  General  Identification. 

Columns  Variable  Description 

1-80  TEXT  Any  identifying  information  for 

the  case 


CARD  2:  Print  Control. 


Columns  Variable 

1-5  IPRI 

6-10  IPRLEV 


Value  Description 

1  All  the  input  data  is 

printed  in  the  same  order 
as  it  is  read 

0  Input  data  print-off 

0  Basic  print  level  (see 

Appendix  C) 

1  Panel  corner  points  are 

printed 


Data , 
Figure 


Format 

8A10 


Format 

215 
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CARD  3:  Mode  Card. 


Columns 

Variable 

Value 

Description 

Format 

1-5 

NPNMAX 

Upper  limit  on  number  of 
panels  the  user  expects 
to  be  generated  (default 
1000) 

215 

6-10 

MSTOP 

0 

Complete  run 

2 

Calculations  stop  after 
SURPAN  for  examination  of 
panel  corner  points  (set 
IPRLEV=1 ) 

CARD  4: 

Number  of  Aerodynamic  Table  Stations. 

Columns 

Variable 

Value 

Description 

Format 

1-5 

NDTAB 

0,2-20 

Number  of  2-D  sectional 
data  table  input  stations 

15 

CARD  5: 

General  Rotor 

Information . 

Columns 

Variable 

Value 

Description 

Format 

1-5 

NB 

1-8 

Number  of  blades 

15 , F10 . 0 

6-15 

RADIUS 

Rotor  radius  as  defined 
in  the  geometry  input 

Note:  RADIUS  must  equal  YFLAP  plus  STY  of  the  last  geometric 

defining  section;  therefore,  RADIUS  sets  the  length  units 
for  geometry  input. 
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CARD  6: 
Columns 
1-10 
11-20 
21-30 


Velocity  Requests . 


Variable 

VCLIMB 

Rotor 

VTIP 

Rotor 

VSOUND 

Speed 

Description 
climb  speed 
rotation  tip  speed 
of  sound 


Format 

3F10.0 


Note:  Any  consistent  units  of  length/sec.  or  dimensionless. 


CARD  7:  Performance  Card. 

Column  Variable  Value  Description  Format 

1-10  CTO  Requested  thrust  coeffi-  2F10.0 

cient?  i.e.r  used  to  set 
the  prescribed  wake  con¬ 
stants 

11-20  RGMX  Estimated  ratial  position  of 

maximum  circulation.  Same 
units  as  RADIUS 


T‘TAKE  DATA 


CARD  8:  Basic  Constants. 


Column 

Variable 

Value 

Description  Format 

1-10 

RFAR 

o 

H 

1 

O 

Ultimate  contracted  radius,  2F10.0 

"A"  in  Landgrebe's  nota¬ 
tion.  Default  value  is  0.78 

11-20 

PSIO 

<^b 

Azimuth  for  the  start  of  axial 
descent  of  the  inner  sheet  at 
the  axis  of  rotation.  In  the 
original  equations,  PSIO=90 
(degrees) 
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CARD  9:  Wake  Azimuthal  Boundaries 


Column 

Variable 

Value 

Description 

Format 

1-10 

PSIFl 

>0 

Azimuth  of  start  of  inter¬ 
mediate  wake  (degrees) 

3F10.0 

11-20 

PSIF2 

>PSIFl 

Azimuth  of  start  of  far 
wake  (degrees) 

21-30 

PSIMRG 

>DPSIT1 

Azimuth  at  which  tip  vortex 
merger  occurs  (degrees) 

CARD  10: 

Wake  Azimuthal  Increments. 

Column 

Variable 

Value 

Description 

Format 

1-10 

DPSIT1 

arbitrary 

Azimuthal  increment  of  tip 
vortex  to  first  blade  pas¬ 
sage  (degrees) 

4F10.0 

11-20 

DPSIT2 

arbitrary 

Azimuthal  increment  of  tip 
vortex  after  first  blade 
passage  (degrees) 

21-30 

DPSIS1 

arbitrary 

Azimuthal  increment  of  inner 
sheet  to  first  blade  pas¬ 
sage  (degrees) 

31-40 

DPSIS2 

arbitrary 

Azimuthal  increment  of  inner 

sheet  after  first  blade  pas¬ 
sage  (degrees) 


Note ;  Unless  NRELAX>0,  all  can  be  set  to  30  degrees  for  pre¬ 
liminary  work. 


CARD  11:  Tip  Vortex  Constants. 


Columns 

Variable 

Value 

Description 

1-10 

TKl 

<0.0 

Pitch  constant  of  the  tip 
vortex  to  first  blade  pas¬ 
sage.  Normally  identified 
as  ki 

11-20 

TK2 

<0.0 

Pitch  constant  of  the  tip 
vortex  after  first  blade 
passage.  Normally  identi¬ 
fied  as  k.2 

21-30 

TLAMDA 

>0.0 

Wake  contraction  constant. 
Normally  identified  as  X 

CARD  12: 

Inner  Sheet 

Constants . 

Columns 

Variable 

Value 

Description 

1-10 

SKTl 

<0.0 

Pitch  constant  of  inner 
sheet  at  r  =  1 . 0  to  PSIO. 
Normally  identified  as 

ki?=i 

11-20 

SKT2 

<0.0 

Pitch  constant  of  inner 
sheet  at  r  =  1 . 0  after 

PSIO.  Normally  identified 
as  k2?=1 

21-30 

SKRl 

<0 . 0 

Pitch  constant  of  inner 
sheet  at  r  =  0  to  PSIO. 
Normally  identified  as  kip_ 
and  taken  as  0.0 

31-40 

SKR2 

<0.0 

Pitch  constant  of  inner 

sheet  at  r  =  0  after  PSIO. 
Normally  identified  as  k^— 


Format 

3F10.0 


Format 
4F10  •  0 


CARD  13:  Tip  Vortex  Formation 


Columns  Variable  Value 


Description 


Format 


1-5  NSEP  M)  The  number  of  tip  panels  I5,2F10.0 

from  the  trailing  edge 
for  which  vortex  shedding 
occurs 


6-15  DELZ 


16-25  DELR 


The  distance  above  the  trail¬ 
ing  edge  that  the  tip  vortex 
passes  through.  Same  units 
as  RADIUS  (Figure  19) 

The  distance  that  the  tip 
vortex  is  shifted  inboard 
from  the  rotor  tip.  Same 
units  as  RADIUS  (figure  19) 


CARD  14:  Wake  Relaxation. 
Columns  Variable  Value 

1-5  N RELAX  >0 


Description 

Number  of  wake  relaxation 
iterations 


Format 

15 


ROTOR  BLADE  ORIENTATION 
(B.C.S.  to  H.C.S.  Transformation) 


CARD  15:  Blade  Orientation. 


Columns 

Variable 

Value  Description  Format 

1-10 

YFLAP 

>0.0  Rotor  flapping  axis  radius  3F10.0 

location.  Same  units  as 

RADIUS 

11-20 

THTA75 

Blade  collective  (degrees) 

21-30 

CONE 

Blade  coning  angle  (degrees) 

Note:  (1)  RADIUS  =  YFLAP  +  STY  of  last  section; 

(2)  see  Section  3,  Figure  3 
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PATCH  GEOMETRY 


Repeat  the  followinq  qroup  of  cards,  16  throuqh  21  (where 


applicable) 

,  for  each 

patch . 

CARP  16:  Patch  Card. 

Co 1 umns 

Vari able 

Value 

Description  Format 

1-5 

I  DENT 

1 

Gives  winq-type  printout  215 

of  pressure  data  (see 

Appendix  C) 

2 

Simplified  pressure  print¬ 
out  (tip  closure  patch) 

6-10 

MAKE 

0 

Reqular  patch  input.  Re¬ 
quires  CARD  17  and  CARDS  18, 

19  and  20  where  applicable 

M 

Automatic  patch  closing  end 
of  Patch  M.  Requires  CARD  21 
only.  See  Section  3.3.2 

CARD  17: 

Section  Card. 

Co  1 umns 

Vari abl e 

Value  Description 

Format 

1-10 

STX 

Location  of  the  section 

6F1 0. 0  , 

n-20 

STY 

loadinq-edqe  in  the  blade 

1  OX, 

21-30 

ST  Z 

reference  frame  (see  Fiqures  415 

6  and  7) 

31-40 

SCAPE 

"  0  Scalinq  factor  on  section 

chord  lenqth 

41-50 

AT,F 

Section  twist  anqle  in^ 

deqrees  I 

:  sec 

Fiq .  6 

51-60 

THETA 

Section  rotation  angle* 

in  deqrees 
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Columns 


Variable 


Value 


Description 


Format 


61-65  INMODE 

0 

-M 

1 

2 

3 

4 

5 


Type  of  section  input 
(see  Figure  8) 

Copies  the  previously  defined 
section  as  originally  speci¬ 
fied;  i.e.,  the  basic  co¬ 
ordinates  before  scaling,  etc. 

Copies  the  basic  coordi¬ 
nates  of  section  M  (absolute 
subscript) 

Input  y,z  coordinates  of 
SECTION.  Requires  CARD  SET 
18A 

Input  x,z  coordinates  of  SECTION. 
Requires  CARD  SET  18A 

Input  x,y  coordinates  of  SECTION. 
Requires  CARD  SET  18C 

Input  x,y,z  coordinates  of  SEC¬ 
TION.  Requires  CARD  SET  18D 

Generates  coordinates  on  an  NACA 
4-digit  section.  Requires 
CARD  20  only 


66-70  NODS  0 

1 

2 

3 

5 


First  or  interior  section  on 
a  catch 

End  of  a  spanwise  region  within 
a  patch  with  continuous  (1)  or 
discontinuous  (2)  slope  on  the 
spanwise  generators  onto  the 
next  spanwise  region  (Figure  11) 

This  section  completes  the 
present  PATCH 

This  section  completes  the 
last  patch  in  the  configuration 


4 
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Columns  Variable 

71-75  NPS 


76-80  INTS 


Value  Description  Format 

Only  active  if  NODS  \  0 

0  Manual  panelling  in  spanwise 

region  just  completed;  defined 
sections  are  used  as  panel 
edges  in  this  region 

^0  Number  of  panels  to  be 

generated  in  the  spanwise 
region  just  completed 

(Only  active  if  NODS  \  0  and 
NPS  0) 

0  j  Form  of  spanwise  interval  spac- 

1  (  ing  for  the  generated  panels. 

2  |  See  Section  3  and  Figure  10 


CARD  SET  18:  SECTION  Coordinates.  (Present  if  INMODE  =  1,  2,  3 
or  4).  In  each  set  use  one  card  per  point.  Insert  NODE 
CARD (S )  19  to  control  panellino  and  to  complete  a  set; 

see  Figures  8  and  9. 


CARD  SET  18A  (INMOOE  =  I) 


Columns 

Variable 

Value  Description 

Format 

1-10 

BY  ) 

y,z  coordinates  of  a  point 

3F10.0 

11-20 

BZ  > 

on  the  section.  The  x- 

21-30 

AX  ) 

stations  are  essentially 
constant  (0.0);  however, 
local  deviations  in  x  can 
be  placed  in  Ax 

- » - 


CARD  SET  18B  (INMODE  =  2) 


Columns 

Variable 

Value  Description 

Format 

1-10 

BX  ) 

x,z  coordinates  of  a  point 

3F10.0 

11-20 

BZ  \ 

on  the  section.  The  y- 

21-30 

AY  J 

stations  are  essentially 
constant  (0.0);  however, 
local  derivations  in  y  can 
be  placed  in  Ay.  Most  rotors 
require  this  form 

CARD  SET 

18C  (INMODE  =  3) 

Columns 

Variable  Value 

Description  Format 

1-10 

BX  ) 

x,y  coordinates  of  a  point  3F10.0 

11-20 

XY  ) 

on  the  section.  The  z- 

21-30 

A  Z  ) 

stations  are  essentially 
constant  (0.0) ;  however,  local 
derivations  in  Z  can  be  placed 
in  Az 

CARD  SET 

18D  (INMODE  =  4) 

Columns 

Variable  Value 

Description 

Format 

1-10 

BX  ) 

x,y,z  coordinates  of  a 

3F10. 0 

11-20 

BY  > 

point  on  an  arbitrary 

21-30 

BZ  ) 

skewed  section 

CARD  19:  Chordwise  Node  Card 


Columns 

Variable 

Value 

Description 

Format 

30-35 

NODEC 

n 

Terminates  a  chordwise 

30X , 31 5 

2  f 

region  havinq  continuous 
(1)  or  discontinuous  (2) 
surface  slope  onto  the  next 
chordwise  reqion  (for  ex¬ 
ample,  see  Figure  B-2) 
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Columns 


Variable 


Value 


Description 


Format 


30-35  NODEC  3 

(Continued) 


-4 


36-40  NPC  0 


-0 


41-45  INTC  0  \ 

1  ( 


Siqnifies  a  TERMINAL  NODE 
CARD  placed  after  the  last 
point  on  a  section 

Neqative  values  initiate  a 
copying  routine.  Use  CARD  19A 
to  define  a  string  of  points 
to  be  copied  over  to  form  part 
or  all  of  present  section. 

The  last  copied  point  is  the 
end  of  a  chordwise  region  with 
the  corresponding  action  ac¬ 
cording  to  the  modulus  of  the 
NODEC  value 

The  CARD  1 9A  must  be  followed 
either  with  another  negative 
NODE  CARD  or  with  further 
basic  points 

Manual  panelling  in  the  chord- 
wise  region  just  completed 
(i.e.,  basic  points  correspond 
to  panel  corners) 

Number  of  panels  to  be  generated 
in  the  chordwise  region  just 
completed 

Form  of  chordwise  interval  spac¬ 
ing  for  the  generated  panels 
(see  Figure  10) 


CARD  I9A:  Copy  Card.  Defines  a  strinq  of  basic  points  to  be 
copied  over.  Required  only  if  NODEC  <  0. 


Columns 

Variable 

Value  Description 

Format 

1-5 

NPCH 

Patch  number 

415 

6-10 

NSEC 

Section  number  within  that 
patch  (local  subscript) 

11-15 

IB  \ 

First  and  last  basic  points 

LB  I 

on  that  section  (local  sub¬ 
script) 
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CARD  20: 


NACA  4-Digit  Section  (INMODE=5)  . 


Columns 

Variable 

Description 

Format 

1-10 

TC 

Thickness/chord  ratio  of  NACA 

F10.0, 

4-Digit  section 

15 

11-15 

NPC 

Number  of  panels  to  be  genera¬ 
ted  across  airfoil  chord  (i.e., 
total  number  of  panels  for 
upper  plus  lower  surfaces  =  2 
x  NPC) 

Note :  A  section  is  generated  having  a  unit  chord  length  and 

having  the  leading  edge  at  0.0,  0.0.  Set  the  STX,STY, 
STZ  and  SCALE  values  on  CARD  17  to  the  required  loca¬ 
tion  and  size  for  the  section. 


CARD  21: 

Automatic 

Patch  Card  (only  if  MAKE^O  on  CARD  16) 

See 

^iaure  12. 

Columns 

Variable 

Value 

Description 

Format 

36-40 

NPC 

Number  of  panels  "chord- 
wise"  on  patch 

35X,  315, 
15X,  315 

41-45 

INTC 

0 

1 1 

V 

Form  of  panel  spacing 
"chordwise" .  See  Figure 

10 

46-50 

KURV 

0 

Flat  tip 

l 

Semicircular  tip  section 

66-70 

NODS 

5 

Last  patch 
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Columns 


Variable 


Value 


Format 


71-75  NPS  0  Basic  points  of  tip  section 

used  as  panel  points  in 
"spanwise"  direction.  NPS 
should  not  be  set  to  zero 
if  NPOO  on  the  tip  section 
(i.e.  ,  on  CARD  19).  If  CARD 
20  is  used  for  the  tip  section, 
NPS  =  0  is  required 

>0  Number  of  panels  to  be 
generated  "spanwise" 

76-80  INTS  0 

11  Spacing  options,  see  figure  10- 

2  \  Inactive  if  NPS  =  0 

3 


AERODYNAMIC  SECTION  DATA 

These  data  tables  are  input  in  the  standard  format  currently 
used  in  the  Rotorcraft  Flight  Simulation  Program,  C-81. 

Repeat  the  following  cards,  22  through  26  (where  applicable) , 
for  each  NDTAB  stations.  If  NDTAB-0 ,  skip  this  entire  card  set. 

CL  and  CD  tables  are  all  that  are  allowed. 


CARD  22:  Copy  Control  and  Defining  Station. 


Columns  Variable  Value 

1-10  ICOPY  0 


ID 


11-20  YR 


Description  Format 

Complete  aerodynamic  I10,F10.0 

tables  are  read  in  for 
this  section 

Sectional  aerodynamic 
data  is  copied  over  from 
previously  defined  sec¬ 
tion  ID 

The  radius  location  of 
the  table  defining  station. 

Same  units  as  RADIUS  (with 
respect  to  axis  of  rota¬ 
tion) 


Note : 


If  ICOPY>0 ,  the  rest  of  the  card  set  (23  through  26)  is 
omitted  for  this  section. 


CARD  23:  Title  and  Control  Card. 


Columns 

Variable 

Value 

Description 

Format 

1-30 

TITLE 

any 

Alphanumerical  title  of 
sets  of  tables 

6A5, 612 

31-32 

NMACH(l) 

3-18 

Number  of  Mach  number 
entries  in  table 

33-34 

NALPHA(l) 

3-125 

Number  of  angle-of- 
attack  entries  in  table 

35-36 

NMACH (2) 

3-18 

Number  of  Mach  number 
entries  in  Cd  table 

37-38 

NALPHA ( 2 ) 

3-125 

Number  of  angle-of- 
attack  entries  in  table 

Note :  Card  sets  24  through  26  are  repeated  as  a  group  two 

times  for  and  C^;  that  is,  in  the  following  descrip¬ 
tions 

K  =  1 - C£ 

K  =  2 - Cd 
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CARD  24:  Mach  Number  Entries. 


Columns 

Variable 

Value 

Description 

Format 

8-14 

MACH (1) 

arbitrary 

Lowest  Mach  number 

7X, 9F10 . 0 

15-21 

MACH ( 2 ) 

m 

Next  highest  Mach  number 

22-28 

• 

MACH ( 3 ) 

• 

H 

ti 

• 

64-70 

• 

MACH (9) 

M 

H 

Note :  Additional  card  may  be  required  with  same  format  to  input 

NMACH(K)  values  of  Mach  numbers. 


CARD  25:  Angle-of-Attack/Coef f icient  Data. 


Columns 

Variable 

Value 

Description 

Format 

1-7 

ALPHAI (K) 

arbitrary 

Angle  of  attack,  degrees 

10F7.0 

8-14 

COI (K, 1 ) 

arbitrary 

Coefficient  at 

MACH (1) 

15-21 

COI (K, 2) 

n 

Coefficient  at 

MACH (2) 

22-28 

• 

COI (K, 3) 

• 

n 

Coefficient  at 

MACH (3) 

• 

64-70 

• 

COI (K, 9) 

ii 

Coefficient  at 

MACH (9) 
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CARD  26: 
Columns 
8-14 

64-70 

Note:  1 

2 


Continued  Coefficient  Data. 

Variable  Value  Description  Format 

COI (K , 10 )  arbitrary  Coefficient  at  MACH(10)  7X,9F7.0 

♦ 

COI (K, 18)  "  Coefficient  at  MACH (18) 

CARD  26  included  only  if  NMACH(K)>9. 

CARDS  25  and  26  repeated  N ALP HA (K)  times  for  each 
angle  of  attack. 
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CARD  1:  TEXT  | 

> 

f 

CARD  2 

. > 

CARD  2 

!:  IPRI,  IPRLEV  j 

l . . 

1:  NPNMAX,  MSTOP  | 

> 

l _ 

CARD  4 :  NDTAB 

\ 

/ 

CARD  5:  NB ,  RADIUS  | 

> 

f 

CARD  6:  VCLIMB,  VTIP,  VSOUND 

. 1 

( 

CARD  7:  CTO,  REMX 

> 

< 

^BASIC  DATA  INPUT  COMPLETE^) 

(a)  Basic  Data. 

Figure  B-l.  Flow  Chart  for  Input  Card  Deck. 
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FROM  BASIC  DATA 


:a) 


If 


CARD  8:  RFAR,  PSIO 


l! 


CARD  9:  PSIFl,  PSIF2 ,  PSIMRG 


CARD  10:  DPSIT1 ,  DPSIT2 ,  DPSIS1,  DPSIS2 


CARD  11:  TKl ,  TK2 ,  TLAMDA ! 


t 


CARD  12:  SKT1,  SKT2,  SKR1,  SKR2 


t 


CARD  13:  NSEP,  DELZ,  DELR 


CARD  14:  N RE LAX 


C 


WAKE  DATA  INPUT  COMPLETE 


0 


(b)  Wake  Data. 
Figure  B-l.  Continued. 
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DEFINE  ANOTHER  PATCH 


I 

i 


i 


(^FROM  WAKE  DATA^ 


CARD  15: 


YFLAP,  THTA75,  CONE 

PATCH  GEOMETRY 


BLADE  ORIENTATION 


>-j  CARD  16:  IDENT,  MAKE 

L  . 


(c)  Blade  Geometry. 


Figure  .  Continued. 
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CARD  22:  ICOPY ,  YR  ! 


>  CARDS  23,  24  ,  25,  26 
CL  AND  CD  DATA 
TABLES 

I  C-81  FORMAT 


<  NDTAB 


=  NDTAB- 


(d)  Aerodynamic  Tables. 


Figure  B-l.  Concluded. 
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SECTION 


© 


Order  of  Input 


Upper  Surface 


Figure  B- 2 . 


Card  Set  18/19.  Section 
for  a  Simple  Airfoil 


Coordinates 


APPENDIX  C:  OUTPUT  DESCRIPTION 


Since  ROTAIR  is  a  preliminary  code,  some  of  the  basic 
printed  output  is  included  until  the  program  has  been  thoroughly 
validated.  Consequently,  only  the  basic  printed  output  is  des¬ 
cribed  here.  Also,  for  the  present  purposes,  the  prescribed 
wake  geometry  is  always  printed  and,  if  NRELAX>0 ,  the  relaxed 
wake  geometry  is  printed  for  each  iteration. 

The  program  output  starts  with  the  program  header  followed 
by  the  user's  TEXT  identification.  If  IPRI  has  been  set  to  1, 
all  the  input  data  is  printed  in  the  order  of  input  and  each 
variable  is  identified  in  parentheses.  The  input  lists  for  the 
Basic  Data  and  Patch  Geometry  are  printed  first.  The  input 
print  is  then  interrupted  by  patch  basic  data  from  subroutine 


GEOMIN . 

This  includes: 

M 

IDENT  KLASS 

KOMP 

NROW 

NCOL 

IPAN 

LPAN 

where  M  is  the  patch  number;  IDENT  is  the  patch  identifier  (1  or 
2) ;  KLASS  and  KOMP  are  not  used  in  the  present  application?  NROW 
and  NCOL  are  the  number  of  panel  rows  and  columns,  respectively, 
on  each  patch,  and  IPAN,  LPAN  are  the  first  and  last  panel  sub¬ 
scripts  on  each  patch. 

If  IPRLEV^l,  the  above  is  followed  by  the  panel  corner 
point  coordinates : 

K  X,  Y.  Z.  X0  Z,  X,,  Z,  X.  Y,  Z  „ 

111  222  333  444 


where  K  is  the  panel  subscript,  and  Xj_,  Y^,  Zi  are  the  normali¬ 
zed  coordinates  in  the  G.C.S..  The  points  are  assembled  in  a 
set  for  each  patch. 

If  IPRI=1 ,  the  input  list  continues  with  a  printout  of  the 
basic  wake  data.  The  prescribed  wake  geometry  is  printed  in  the 
same  form  as  for  the  HOVER  code.  However,  actual  normalized  co¬ 
ordinates  are  defined  with  respect  to  the  G.C.S. 
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The  printout  of  the  pressure  distribution,  etc.,  is  as 
follows: 

(a)  MODE  =  1.  The  pressure  distribution  is  printed  according 
to  the  type  of  patch  (IDENT) .  If  IDENT  -  1,  the  following  is 
printed  for  each  column  of  panels. 

K  X  Y  Z  DUB  VX  VY  VZ  V  C  X/C  C* 

jj  ir 


where  K  is  the  panel  subscript;  X,  Y,  Z  is  the  control  point 
location,  and  DUB  the  doublet  value;  VX,  VY ,  VZ  and  V  are  the 
velocity  vector  components  and  magnitude;  Cp  is  the  usual  pres¬ 
sure  coefficient,  X/C  is  the  location  of  the  point  relative  to 
the  local  chord  line,  and  C*  is  the  pressure  coefficient  normali¬ 
zed  by  VTIP. 

Each  column  is  followed  by  the  integrated  section  force  and 
moment  coefficients  defined  in  the  G.C.S. 

CIRC  CX  CY  CTS  CTSGAM  CLGAM  CL  CD  CMX  CMY  CQS 


where  CIRC  is  the  normalized  section  circulation;  CX,  CY  are 
the  force  coefficients  along  the  X-  and  Y-axes  of  the  G.C.S. ;  CTS, 
CTSGAM  are  the  section  thrust  coefficients  by  pressure  integra¬ 
tion  and  circulation,  respectively;  CLGAM,  CL  are  the  section 
lift  coefficients  by  circulation  and  pressure  integration,  res¬ 
pectively;  CD  is  the  section  drag  coefficient  by  pressure  inte¬ 
gration;  CMX,  CMY  are  the  moment  coefficients  about  the  X-  and 
Y-axes  of  the  G.C.S.;  and  CQS  is  the  section  torque  coefficient. 
CHORD  is  the  local  chord  length  with  leading  edge  at  XLE ,  YLE, 
ZLE.  The  coordinates  are  given  in  the  same  units  as  originally 
specified . 

For  type  2  patches  (IDENT  =  2) ,  the  pressure  printout  is 
shortened  by  omitting  X/C,  C*. 

A  summary  of  the  patch  force  and  moment  coefficients  nor¬ 
malized  by  the  reference  quantities  in  the  basic  data  is  printed 
after  each  patch  and  the  total  force  and  moment  coefficients  are 
printed  after  the  last  patch. 

If  NDTAB> 2 ,  then  the  profile  losses  calculated  by  the  user- 
supplied  two-3imensional  tables  are  printed  in  the  following 
manner. 


i 
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COL  Y  MACH  CL2D  CDO  CXO  CTSO  CQSO 

•  ••  •  •  •  •  • 

•  •  *  •  •••  • 

•  •  •  •  •••  • 

COL  is  the  column  identified;  Y  is  the  normalized  radial 
position  of  the  section  in  the  G.C.S.;  MACH  is  the  section  Mach 
number;  CL2D  is  the  interpolated  table,  CL  (should  equal  the 
surface  singularity,  CL);  CDO  is  the  corresponding  interpolated 
profile  drag  coefficient;  CXO  is  the  resultant  sectional  inplane 
profile  force  coefficient;  CTSO  is  the  resultant  sectional  pro¬ 
file  thrust  coefficient  loss;  and  CQSO  is  the  resultant  sectional 
profile  torque  coefficient. 

Finally,  the  integrated  total  rotor  performance  properties 
(based  on  pressure  integration)  including  the  figure  of  merit 
are  printed. 
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APPRENDIX  D:  SAMPLE  INPUT  DATA  CASES 


CASE  Is  HYPOTHETICAL  ROTOR  GEOMETRY 


Geometric  Description  (scaled  by  R) : 


Tip  Taper  Ratio 

= 

1.5 

%  Radius:  Start 
Taper 

Of 

= 

85 

Blade  Chord 

Root  to  Start  of 
Taper 

_ 

.  06 

Tip  Chord 

- 

.04 

Sweepback  of  c/4 

= 

0 

O 

CM 

%  Radius:  Start 
Sweep 

o  f 

= 

85 

Anhedral  of  c/4 

= 

10° 

%  Radius:  Start 
Anhedral 

of 

= 

85 

Airfoil 

= 

NACA  0012 

Linear  Twist  about 
c/4 

_ 

-10° 
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TABLE  COPY  CARO 


APPENDIX  E 


SAMPLE  OUTPUT  DATA  CASE 


ANALYTICAL  Hfc I  MOOS  *  I NC • 

PROGRAM  R[)TA]R 

SUklACf  POT  t  M  I  Al  Kl)n  ANALYSIS  M  t  T  HOD 
f  OR  ROTORS  IN  HOVtR  Ok  CL  MB 

UHSttA  ROT  Ok 

$A*ev4«t««*$«***$e$t**$$«*t$**t»$$*£«**$*$*t«**et* *«*$$$* 
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Cf'i-  v  ClSl,  COO  CXU  CTSJ  COSJ 


LIST  OF  SYMBOLS 


B,C 


c 

C 

P 

CL 

Cm 

E 

i'i'k 

N 

n 

R 

S 

s 

V 

v 

W 


x,y,z 


wc 

w 

climb 


u 


Velocity  potential  influence  coefficients  for  constant 
source  and  doublet  distributions,  respectively,  on 
quadrilateral  panels,  Eq  (5) 

Airfoil  chord 

Pressure  coefficient 

Lift  coefficient 

Moment  coefficient 

Quantity  in  Eq  (5) 

Orthogonal  unit  vector  system  defining  the  axes  of  the 
Cartesian  coordinate  system  fixed  relative  to  the  blade. 
Figure  2 

Number  of  panels 

Surface  unit  normal  vector  directed  into  the  flow  field 

Position  vector  of  a  point  relative  to  a  point  on  the 
axis  of  rotation 

Blade  surface 

Surface  distance 

Velocity 

Perturbation  velocity 
Wake  surface 

Cartesian  coordinates  in  th^  blade-fixed  frame 

Normalized  climb  speed 

Climb  speed 

Body  rotation  vector 

Finite  increment 

Doublet  strength 
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Total  velocity  potential 
Perturbation  velocity  potential 
Source  strength 


$ 

<P 
o 

aR'-M  Source  components  due  to  unit  pitch  rotation  and  unit 
translation,  respectively,  Eq  (7) 

V  Gradient  operator 


Subscripts 

J,K  Values  on  panels  J,  K,  respectively 

U , L  Upper,  lower 

W  Wake 

00  Reference  onset  condition 
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